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, I explicitly calculate the anomalous dimensions and splitting functions governing the evolu- 

. tion of the parton densities and structure functions which result from the running coupling BFKL 



equation at LO, i.e. I perform a resummation in powers of ln(l/x) and in powers of (3q simultane- 
ously. This is extended as far as possible to NLO. These are expressed in an exact, perturbatively 
"q^' calculable analytic form, up to small power-suppressed contributions which may also be modelled 
Q^l to very good accuracy by analytic expressions. Infrared renormalons, while in principle present in 
r~| ! a solution in terms of powers in a<j((5^), are ultimately avoided. The few higher twist contribu- 
tions which are directly calculable are extremely small. The splitting functions are very different 
^ ■ from those obtained from the fixed coupling equation, with weaker power-like growth ~ x~^ '^^, 
■ which does not set in until extremely small x indeed. The NLO BFKL corrections to the splitting 
functions are moderate, both for the form of the asymptotic power-like behaviour and more impor- 
tantly for the range of x relevant for collider physics. Hence, a stable perturbative expansion and 
predictive power at small x are obtained. 
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1. Introduction. 

Small X physics has been a particularly active area of particle physics research in the past few 
years, driven largely by the first data for x < 0.005 being obtained by the HERA experiments [1] 
[2]. However, as well as the need to describe this HERA data correctly, it will also be extremely 
important to understand the correct way of calculating physics at small x in order to interpret the 
results coming from the LHC in a truly quantitative manner. For example, for the production of a 
particle of mass ^ lOOGeV the typical value of x probed (at central rapidity) is 0.005, but values 
up to two orders of magnitude in either direction will also have an almost equally large influence.^ 

The potential complication at small x is that the splitting functions and coefficient functions 
governing the evolution of parton distributions and their conversion to physical quantities have 
terms in their perturbative expansions which behave like a" ln™'(l/x), where m can reach up to 
n — 1. Therefore, as the power of the coupling increases the powers of = ln(l/x) also increase, 
and rapid perturbative convergence is not really guaranteed if ^ ^ l/a^ i.e. ~ 5. This problem 
is not really diminished at the LHC, where the coupling is likely to be smaller than at HERA, 
since the parton distributions to be used will be those measured at HERA at much lower scales 
and evolved up to LHC scales. This question of large ln(l/x) terms is in principle addressed by 
the BFKL equation [4], which is an integral equation for the unintegrated 4-point gluon Green's 
function in the high energy limit. This sums the leading high-energy, or in the DIS case, small-x 
behaviour, which is dominated by the gluon, and thus allows the extraction of leading ln(l/a;) terms 
for relevant quantities, such as splitting functions. 

Hence, a major point of debate of the past decade has been whether the standard DGLAP 
approach based on renormalization group equations and conventionally ordered simply in powers 
of as{Q^), or the BFKL equation, which sums leading logarithms in is most effective way 

of dealing with small x physics (most particularly structure functions), and/or whether the two 
approaches need to be combined in some way, and if so, how? While the conventional DGLAP 
approach has been relatively successful, it does have some significant problems (which are often 
overlooked): a valence-like, or even negative input gluon leading to a strange low-Q^ Fl{x,Q'^); 
undershooting of the data systematically for x ~ 0.01 at the highest when a global fit is 
performed; and apparent instability at small x order-by-order in ag up to NLLO [5].^ 

Nevertheless, the BFKL equation did not seem to help these problems. The original LO 
BFKL prediction of a behaviour of the form x~^ for structure functions and splitting functions 
at small x, with A ~ 0.5, was clearly ruled out long ago. A combination of the two approaches, 
using the BFKL equation to supplement the Altarelli-Parisi splitting functions with higher terms 

^ For an illustration of the x and of parton distributions sampled at the LHC see fig. 1 of [3] . 
^ Of course the full NNLO splitting functions are not known, but good estimates are available [6] based 
on calculation of moments in [7] . 
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of the form a"+^ ln"'(l/a;) was originally successful (so long as one avoided factorization scheme 
ambiguities by working in physical quantities) [8], but this success is not possible to sustain with 
the most recent data [9] [10]. Moreover, the subject was thrown into confusion by the calculation 
of the NLO correction to the BFKL equation [11] [12]. The results of this calculation were not very 
encouraging. Ignoring the running of the coupling at NLO, i.e. proceeding with the same sort of 
calculations as at LO but including the scale-independent NLO correction to the kernel, one obtains 
the "intercept" for the splitting function power-like behaviour, shifted from A = Aln2as to 
A = 41n2as(l — 6.50;^). This is clearly a huge correction, and implies the breakdown of the 
perturbative expansion for this quantity. More serious than this intercept is the power series for 
the splitting function, which may be calculated even taking into account the renormalization and 
scale dependence introduced at NLO. Expanding this out formally to NLO in ln(l/x) one finds that 
it is is dominated by the NLO corrections at all values of x below about x = 0.01. For example, 
using the formulae in [11] the first few terms in the power series for P{x) go like 

xP{x, Q^) =as + 2Aa1i^/Q + 2.1a^^Vl20 + • • • 

- as(0.43a^ + 1.6a^^ + ll.l ali^ / 2 + IS.Sa^^V^ + 39.7a^^^/24 + 169.4a^^Vl20 + •••), 

(1.1) 

where £^ = ln(l/x) and ctg = agiC^). Clearly, the size of the coefficients more than compensates 
for the extra power of ^^(Q^), particularly at low where the perturbative analysis of structure 
function evolution often takes place. 

Hence, this NLO correction left open the whole question of how to address the evolution of 
structure functions at small x. There has been considerable progress on the stability of the solutions 
to the BFKL equation in the intervening time. One major development was the observation that 
the resummation of double logarithmic terms in the transverse momentum is necessary in order 
to eliminate collinear divergences. This renders the intercept of the BFKL equation stable [13], even 
when ignoring the renormalization scale dependence. This initial idea has been further developed 
in [14] [15] [16], where the effect of running coupling is also considered in these later papers. This 
development is particularly important for the case of so-called "single scale" processes where both 
ends of the gluon Green's function are at high scales (not necessarily the same) where without this 
collinear resummation, all calculations are badly behaved over the full range of energy, not just in 
the asymptotic limit. 

However, for the type of situation embodied by DIS, where one end of the gluon Green's 
function is at some low non-perturbative scale, the factorization theorem simplifies the problem. 
Although the growth of the coupling at low scales actually renders the solution of the BFKL 
equation formally divergent when the renormalization of the coupling is encountered, as realized 
as long ago as [17] and studied in detail in [19], all the uncertainty and indeed all the effects of the 
low region are absorbed into the overall normalization of the gluon, leaving the evolution and 
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coefficient functions for hard scattering cross-sections calculable. However, these perturbatively 
calculable quantities are affected by the running of the coupling, and it was argued in [20] that 
the effective result was as if the usual LO BFKL splitting functions should be evaluated at an 
x-dependent scale, which grows with decreasing x, due to increasing diffusion into the ultraviolet, 
leading to a decrease in the coupling. Hence, the effect of running coupling totally transforms 
the more simplistic LO BFKL results, making overall normalization of quantities incalculable, 
but moderating the effect of those governing the evolution in Q^. This moderation of the LO 
quantities also translated into a moderation of the effects of NLO corrections, leading to a much 
improved stability of the perturbative expansion, even without recourse to the type of resummation 
in [13]- [15]. Indeed, for this case of deep inelastic scattering further resummation of this type is 
redundant. These modified BFKL contributions to the splitting functions, when combined with 
the conventional LO-in-ag contributions, also led to improved fits compared to the usual DGLAP 
approach [20] and a more sensible prediction for Fl{x,Q^). This concept was put on a firmer 
footing in [21] where an expHcit calculation of the BFKL splitting functions in powers of Poas{Q^), 
i.e. a resummation of running coupling contributions, was outlined, and it was seen that over a 
wide range of the x — range (including the HERA range) the previous hypothesis was largely 
correct, and precise results were also obtained outside this range. 

The purpose of this paper is to explain in detail, and expand upon the results of this previous 
letter, i.e. to present in full the calculation of splitting functions and coefficient functions for 
deep inelastic scattering obtained from the BFKL equation (both LO and NLO) and incorporating 
running coupling contributions to all orders. Explicitly, while the usual BFKL equation presents 
an expression for these quantities which sums the leading power of ^ at each power in a^, I will 
extend this by producing expressions which also include the leading power of /3o at each power of 
as{Q'^) and e.g. 

OO n — 1 

xPgg{x,Q^) = Y.Y1 «™<(Q')r-'-'"/3o"^, (1-2) 

n=l m=0 

though the formal divergence of the series will complicate this form a little. This presentation will 
begin, in section 2, a brief review of the standard solution to the BFKL equation at LO, and then 
a detailed presentation of the solution at LO with running coupling. This will result in a solution 
for the gluon splitting function in an analytic form up to a small, unambiguous, correction of the 
form A^/Q^ (which is not higher twist) which may be modelled by an analytic function to excellent 
accuracy. Despite the integration over the infrared region when solving the running coupling BFKL 
equation there is no ambiguity in this splitting function. Next, in section 3, will follow a discussion 
of some possible higher twist contributions at small x. It is argued that these are may be much 
smaller than generally supposed, though the possibility of some large power-suppressed corrections 
(not necessarily higher twist) is left open. In section 4 1 discuss the solution of the BFKL equation at 
NLO, defining precisely what I mean by the "NLO BFKL splitting function", and showing that the 
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NLO corrections for the gluon splitting function are moderate. In section 5 I consider real physical 
quantities, i.e the structure functions. Firstly, I calculate the quark-gluon splitting function and 
coefficient functions, and then consider the rather more direct physical splitting functions [22]. I 
also consider how far one can calculate to NLO, defining a "nearly NLO" physical splitting function 
Pll{x,Q^)- The stability of the perturbative expansion is examined in detail, and seen to be very 
good. Finally, in section 6 phenomenology is briefly touched upon, and I present a summary and 
my conclusions. 

2. BFKL Equation at LO. 

The BFKL equation for zero momentum transfer is an integral equation for the 4-point, trans- 
verse momentum-dependent gluon Green's function for forward scattering in the high energy limit, 
f{ki, k2,as, N), where is the McUin conjugate variable to energy. In the case of DIS the second 
momentum k2 is put equal to some non-perturbative scale Qq, we let fei = k, and N becomes 
conjugate to x. In order to obtain a structure function we attach the non-perturbative bare gluon 
distribution gB{N,Ql) to the non-perturbative end of the gluon Green's function and convolute a 
hard scattering cross-section h{Q^ /k'^ ,as, N) to the perturbative end. 

In this section I will illustrate the effect that introducing the running coupling into the BFKL 
equation has. In order to do this I will first begin with a brief presentation of the fairly simple tra- 
ditional case of fixed coupling before moving to the far more complicated case of running coupling. 
As will be seen, the introduction of renormalization, and hence running of the coupling, which is 
necessary except in the artificial model of no consideration beyond LO, completely changes not 
only the detail of the information one is able to extract from the BFKL equation, but also what 
type of information one is able to extract. 

2.1. Fixed Coupling. 

We simplify matters by working in moment space, i.e. defining the moment of a structure 
function by 



and similarly for the parton distributions (scaled by x). Doing this the BFKL equation is 



where f{k'^,as/N) is the unintegrated gluon four-point function, fi{k'^,Qo) is the zeroth order 
input, as = (3/7r)as, and the LO kernel is defined by 




(2.1) 




(2.2) 




(2.3) 
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It is convenient to define the input by fi{k'^,Qo) = S{k^ — Qq). In fact in the leading twist 
factorization theorem this is the unique definition, and Qq is really just a regularization which we 
let ultimately. Going beyond this approximation the dependence on Qq tells us about the 
higher twist due to the intrinsic transverse momentum of the gluon, and we will discuss this in 
section 3. The "gluon structure function" is now given by 

g{Q^N) = -^f{N,k^Ql)xgB{N,Ql), (2.4) 

where gsiNjQ^) is the bare gluon density in the proton which implicitly absorbs the collinear 
divergences in /(A;^). The BFKL equation is most easily solved by taking the Mellin transformation 
to 7-space, i.e. 

/(7,iV)=/ dk\kY'-y{k\N), (2.5) 
Jo 

where it reduces to 

/(7,iV) = fi{l,Ql) + (a./iV)xo(7)/(7,iV), (2.6) 
where /(7, Qo) = 6xp(— 7ln((5Q)) and x{l) is the characteristic function 

Xo(7)=2V'(l)-V'(7)-V'(l-7)- (2.7) 
A little simple manipulation leads to the expression 

This inverse transformation has a leading twist component given by the contribution of the leading 
pole at 1 — {o:s/N)xo{7) = 0, and the solution is 



The anomalous dimension jo{as/N) may be transformed to x-space as a power series in ag ln(l/a;), 
and has a branch point at TV" = A = 41n2Qs (at which 7 ^ |) leading to asymptotic small x 
behaviour for the splitting function 

xP^^ix,a,)^ (^^^)3/2 ^ (2-10) 

In a similar fashion, assuming that the leading small x behaviour is dominated by the perturbative 
physics rather than by gsiQa, TV), one can transform to x-space the normalization /j^^^^^i (^^^ 
finding that this leads to a gluon normalization xg{x) oc g)i/2 • 



2.2. Running Coupling. 

Beyond strict leading order it is impossible to ignore the running of the coupling. At NLO 
ultraviolet regularization is required, resulting in a correction to the LO kernel of the form 
—(5oas{lJ'%)lyi{k^/lJ'%)KQ{q^,k^), where iJ:r is the renormalization scale which must now be in- 
troduced. Hence, it is unrealistic to simply use the LO kernel without considering the influence 
of such a correction. An obvious way in which to incorporate such a term is to simply use the 
running coupling constant evaluated at the scale k"^ in the previous LO BFKL equation. Since this, 
or something similar, is unavoidably forced upon us at NLO, it seems sensible to consider the fixed 
coupling LO BFKL equation as just a model which would apply in a conformally invariant world, 
and more realistically to work with the BFKL equation with running coupling [23] [24] [17] [18] from 
the beginning. Doing this we obtain 

fie,QlaAk')/N) = h{k\Ql) + ^ r ^Ko(«^fe^)/(g^), (2.11) 

Jo 1 

where 

a, = l/(/3oln(fcVA')), (2.12) 

Po = (11 — 2Nf/3)/{4:TT), and Nf is the number of active flavours. 

One can solve this equation in the same type of way as for the fixed coupling case, i.e. take 
the Mellin transformation, but now with respect to (fc^/A^). It is most convenient first to multiply 
through by ln(fe^/A^), in which case one obtains 

where Po = (tt/Jq / 3) . Hence, the inclusion of the running coupling has completely changed the form 
of our double Mellin space equation, turning it into a first order differential equation. This has a 
profound effect on the form of the solutions. The equation may easily, if formally, be solved giving, 

/(7, N) = eM-M7)/iPoN)) H ^M^^J^^Qo) exp(Xo(7)/(^oiV))c/7, (2-14) 

where 

Xo(7) = |^\o(7)d7^ (2V(l)(7-^)-ln(^;^^)). (2.15) 

-^0(7) — ^ 111(7) at 7 = and hence exp(— Xo(7)/(/3oA^)) has a branch point at 7 = 
(exp(— Xo(7)/(/3oA^)) -y-'^/PoN-^ wit\\ similar branch points at all negative integers. It is eas- 
iest to choose each of the cuts along the negative real axis. exp(Xo(7)/(;5oA^)) has similar branch 
points at every positive integer, and it is easiest to choose these cuts along the positive real axis. 
This means that the integral in (2.14) is ambiguous due to the available choice in avoiding the 
cuts. This ambiguity can only really be removed by regulating the Landau pole in the definition of 
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the coupling. However, this introduces model dependence, and also makes analytic progress rather 
more difficult, so I simply accept this ambiguity for this function.^ 

In order to simplify (2.14), and introduce factorization we trivially rewrite it as 



/(7,iV)=exp(-Xo(7)/(^oiV)) 



exp(Xo(7)/(^oA^))d7, (2.16) 



In the region of 7 = the integrand in (2.16) is oc j^/I^oN ^ ^j^g integral of this from ^ 7 is 
oc 7i+i//3o^. Hence, the leading singularity in the 7 plane for exp(— Xo(7)/(;9o-^))) is cancelled 
by the integral from — 7 of this integrand [18], and the new leading singularity is at 7 = — 1. 
Since G{Q^,N) is obtained by an inverse Mellin transformation with respect to Q^/A^, the part 
of (2.16) coming from the integral from to 7 will behave like A^/Q^ (actually Qq/X^ as we will 
see later). Hence, disguarding this power-suppressed correction, which will be considered in some 
detail in section 3, we keep only the first term in (2.16), obtaining for the gluon distribution 



1 

g{Q',N) = — L - exp(7ln(gVA') - Xo{^)/{PoN))dj 



2 — «oo 

■)2/a2\ I V ix,\na ATWJx, ^ irSi 



f 

Jo 



eM-lHQo/^')+Ml)/iPoN))djgB{QiN) 
= gE{Q\N)gi{QlN)gB{QlN). 

(2.17) 

Therefore, we have factorization up to well-defined corrections of 0{Qq/Q^), which genuinely 
do vanish as Qq ^ (see section 3). As mentioned, exp(Xo (7)7(^0-^))) contains singularities at 
all positive integers, and Qi{Q'q-,N) is not properly defined, since the integrand has singularities 
lying along the line of integration. However, since this factor is independent of Q^, it does not 
contribute at all to the evolution of the structure function. It is also divergent as ~^ ^1 
as usual in the factorization theorem these divergences are implicitly cancelled by gB{Qo,N), and 
we can imagine the ambiguity to be cancelled in the same manner. So the overall normalization 
is incalculable, but there is a calculable function Qe{Q^,N) whose form is determined by the 
singularities of exp(— Xo(7)/(;9o-^)) in the 7 plane. This also leads to a fundamental difference 
between the cases of the fixed and running couplings. Whereas previously the leading singularity 
was a pole at {as/N)x{'j) = 1, i.e. at 7 — > | as A?^ ^ 41n2as, now the leading singularity is an 
cut at 7 = 0: there is no power- like behaviour in Q^. Similarly, the branch point in the N plane at 
41n2as has become an essential singularity at AT = 0: there is no power-like behaviour in x in the 



^ The problem due to the Landau pole is illustrated using an alternative method of solution in 
[19]. In this paper the solution of the equation where the NLO coupling effect is left simply as 
—poas{iJ.^ji)ln{k^/lj,^)Ko{q^,K^) rather than resummed is also considered. This does not improve the 
situation, i.e. an ambiguity in the solution remains even in this case. 
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evolution factor for the gluon. The introduction of the running of the couphng has changed the 
character of the solution completely. 

One can now proceed with the solution to the LO BFKL equation by acknowledging that the 
only real information contained in Qe{N,Q^) is on the evolution of the structure function, i.e. 
defining 

d\nG{NM^) dlng^(jV,Q^) 2 

dln(Q2) = dln(Q2) =7..(A^,Q )■ (2-18) 

Qe{N-,Q'^) gives us an entirely perturbative effective anomalous dimension governing the evolution 
of the gluon structure function. The usual technique for solving for Qe{N ^Q"^) is to expand the 
integrand in (2.17), about the saddle-point. This results in a contour of integration parallel to the 
imaginary axis, with real part | for the small x solutions, see fig. 1. Using this results in an 
anomalous dimension 

oo 

^gg{N,Q'')=^o{as{Q'')/N) + Y,{-(ioas{Q'')rin{as{Q'')/N), (2.19) 

n=l 

i.e., the effective anomalous dimension is the naive leading-order result with coupling at scale 
plus an infinite scries of corrections in increasing powers of — /3oas(Q^) [20]. However, each 
of the j{as{Q'^)/N) is singular at = X{Q'^), and the power of the singularity increases with 
increasing n. Hence, although the series for the resulting splitting function is in the small quantity 
as{Q'^)Po, the accompanying coefficients are progressively more singular as a:; ^ 0. The saddle- 
point approximation is therefore not a reliable result as a; ^ and explicit investigation reveals that 
it is only really quantitatively useful when as{Q^) ln(l/a;) is so small that the effective anomalous 
dimension is effectively the LO in ag part, xPgg{x) = as{Q^)[20]. This translates into x <^ 0.01 
in the HERA range. Therefore the calculations of the anomalous dimension which rely on an 
expansion about the saddle-point, i.e. the conventional expansion in decreasing powers of ln(l/x) 
at fixed power of ag, leads to very inaccurate and misleading results for small x. This instability 
is not surprising. If one examines the integrand along the saddle-point contour of integration one 
finds that it is very different from the Gaussian form the saddle-point method assumes [20]. Also 
this is an expansion obtained from approaching 7 = | and in terms of functions of N which are 
singular at N = X{Q^), whereas wc know that the full solution no longer sees these points as 
anything special. In fact, the known singularity structure of the integrand implies that 7 = is 
the point on which to concentrate. 

This suggests an alternative method of solution for the anomalous dimension. In order to 
concentrate on this leading singularity we may move the contour of integration to the left and 
simultaneously use the property that the integrand dies away very quickly at infinity (for 5?e7 < 1 /2) 
to close the contour so that it simply encloses the real axis for 7 < fig. 1. It is then useful to 
express Xo(7) in the form 

00 

Xo(7) = 1/7 + 1; 2C(2n + 1)7'", (2.20) 

n=l 
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which is, however, only strictly valid only for I7I < 1. Doing this we may write 

^0(7) = ln(7) +7^ + 1; 2^|j^7^"+\ (2.21) 

n=l 

and the integrand for Qe{N, Q^) becomes 

^-i/(^oiV)-i3^pj^^,_ _i_(^^^|j^^^2„+i)^^ (2.22) 



where t = ln((5^/A^) and a„ = 2^(2n + l)/(2n + 1). The contribution to the integral from 
— > —00 + ie is now the same as that from —00 — — > up to a phase factor, and we may write 

(2.23) 

where the integral has to be understood as an analytic continuation, since there are singularities 
along the real axis, and strictly speaking the integrand is well defined only for 7 > —1. Since the 
factor of exp(7t) is present this latter point leads, in principle, to an error of order exp(— t), i.e. 
0{h? /Q'^) into the value of QE{N,t). This will be discussed in more detail below. 
In order to evaluate the above integral it is convenient to let y = ^yt, resulting in 

gE{N,t) = -sm(-J^]e^J- 



i^/^^"^) J ^y-^/(^°^)-^exp(y)exp(^-^X;a.(yA)2"+ijd2/. 

The latter exponential may be expanded as a power series in y/t and the each term in the integral 
then precisely evaluated using the standard result that 

(_l)"r(-l/(^oA^) + n)= / y-V(/3oAf)-i exp(y)y"dy. (2.25) 

•J —00 

Hence, we may formally write 



plus an error of 0{A'^/Q^). We note that we could have reached this final expression (2.26) in a 
slightly more rigorous manner. After performing the expansion of ^0(7) in (2.21) we could have 
produced a well-defined integral in (2.23)by taking the lower limit of integration to be —1 + e so 
that the expansion is valid over the region of integration. This would mean that there is region of 
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integration 7 < — 1 absent, which due to the factor of exp(7t) would mean a missing contribution of 
0{h? IQ"^). This new hmit of integration would result in the lower limit of —t in (2.24) and (2.25) 
and consequently we would obtain incomplete gamma functions ^{—\/{^qN) + n,t) rather than 
r(-l/(^oA^) + n). However, 7(-l/(^oiV) +n,t)= T{-l/{PoN) + n) + 0{h?/Q''), so disguarding 
the contributions of 0{A'^/Q^) we regain (2.26), which is formally equivalent to (2.23), but we have 
seen explicitly the origin of the intuitively obvious 0(A^/Q^) corrections to (2.26). 

The result (2.26) was first noted in [25], and was simplified by using the relationship that as 
N ^0, {r{-l/{PoN)+n)/T{-l/{PoN))) (-1/(^0-^))"- However, it is important to notice the 
more general result that for all N 

where 

n-l 

A„(-l/(^oiV)) = ^ (-l)-d„„(^oiV)-"+-, (2.28) 

m=0 

and dmn are positive coefficients and don = 1- Explicitly the first few An{—l/{f3oN)) are 
Ai(-l/(/3oiV)) ' ^ 
A2(-l/(^oA^)) 



(/3oiV), 

1 V r I 



{l3oN)J V(/3oiV) 



1 \" 1 V J 1 



A3(-l/(/3oA^))= -3 +2 



(/3oiV)y V(/5oiV)y V(/3oiV) 



A4(-l/(/3oA^))= 7^ -67^ +11 7^ -6 



(/?oAr)y V(/3oA^)y V(/3oiV)y V(/9oiV) 

A5(-l/(/3oiV)) = ( -lofy^^Tvl +35f75\7vl '^^(l^^ +24 



(/?oAr); V(/?oiv)y V(/?oiv)y VC^iv); 

(2.29) 

These functions oscillate a great deal and only approach the asymptotic values of l/(/3oA^)" at low 
values of N which decrease with increasing n. The comparison of A4(— l/(^o-^)) with l/{PoNY is 
shown in fig. 2, and illustrates this feature clearly. 

Ignoring the common factor of — sin(7r/(/3oAf))r(— l/(/3oA^)) exp(— 7£;/(/3oiV)), which has no t 
dependence, and is irrelevant for the calculation of the anomalous dimension, 

GE{N,t) = t^/CPoN)L + ^^„(i/(^oiV))i-"A„(-l/(^oiV))) (2.30) 

^ n=3 ^ 

where the A„ are simply calculable from the expansion of exp^— l/(^o-^) X^^i o„(y/t)^"'"'"-^^ . 
The common factor of t^/^PoN) jg ^j^g well-known double-leading-log result coming from just the 
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LO as{Q^)/N part of the anomalous dimension. Multiplying this we have an expansion as a power 
series in 1/t or equivalently in as{Q^). In fact 

n-l 

i-"A„(-l/(^oiV)) = (a,(Q2)/iV)- ^ dmni-Poas{Q^)r{a,{Q^)/N)-^. (2.31) 

m=0 

This explicitly demonstrates that we obtain a set of running coupling corrections to a LO result, 
i.e. in solving the BFKL equation we are now obtaining not only the leading power in 1/N 
(corresponding to the leading power of ln(l/a;)) at each order in as{Q^), but we also obtain the 
leading power in /3o at each power of as(Q^) and 1/N. Substituting this type of expansion into 
(2.30), putting the resulting expression for QE{N,t) in (2.18) and expanding in inverse powers of 
one obtains an expression for the anomalous dimension as a power series in asiQ"^), where at each 
order we have the leading divergence in 1/N plus a sum of running coupling correction type terms. 
With a little work one may regain the whole leading 7o(as(Q^)/A^) (though it is necessary to keep 
some subleading terms in the A„ to do this), along with a tower of terms which are subleading 
in powers of /3oas(Q^) to this leading anomalous dimension - one obtains all the corrections to 
this naive LO anomalous dimension due to the running of the coupling, i.e. the whole of (2.19) is 
regained, but ordered in powers of UsiQ"^) rather than in /3oas(Q^). 

The general features of this full, running coupling BFKL gluon Green's function and consequent 
anomalous dimension may be appreciated quite easily. The important fact to note is that although 
the A„(-l/(^oA^))) (l/(/^o^))" as iV ^ 0, the function oscillates a great deal with 1/{(3qN), 
and remains much smaller in magnitude than this asymptotic form until very small N , roughly 
until 1/iV > n. This coupled with the accompanying factor of means that for reasonable i.e 
t ^ 4 - 5 (Q2 ^ IGeV^), only the first 5 or so terms in (2.30) make a significant contribution for 
N > 0.25. Hence, to a very good approximation 

g.(N,t) = (i _ _m_A3(-l/(A^)) - -|i|,A,(-l/(A^))), (2.32) 

and in fact the smallness of the coefficient makes even the term almost negligible in this case. 
GsiNji) initially grows as N falls due to the t^/^l^oN) ^^j.^ However, for N ~ 0.6 the negative 
contribution from the term starts to become significant and ultimately drives the gluon structure 
function to negative values. The result is shown in fig. 3. dQE{N,t)/dt may simply be evaluated 
also using (2.30), and shows the same general shape, but does not become negative until a slightly 
lower value of N as also seen in fig. 3. Hence the anomalous dimension develops a leading pole at 
a finite value of N, given by 

,^2a^(_i 3_ ^\ 

3(;3oiV)V(/3oAr)3 C^^NY {PoN)J- ^ ' ' 
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This result is accurate to better than 10% even at ~ IGeV^, and is much better at higher Q^, 
the right-hand-side receiving corrections formally of 0{l/{t'^^QN^)), but which are numerically 
small. The value of N for this leading pole is shown as a function of t in fig. 4, and for the sort of 
values of t relevant at HERA is ~ 0.25. Going to N < 0.25 higher order terms in (2.30) become 
important, and the positive l/{{PoN)H^)Afi{-l/{PoN)) term absent in (2.30) pulls QE{N,t) back 
to positive values, and another pole, with opposite sign residue, appears in ^gg{N,t). At even lower 
N the analytic expression eventually breaks down, as discussed below, but numerical results show a 
series of poles becoming closer together. Nevertheless, the position of the leading pole is essentially 
determined by the first handful of terms in the power series in UsiQ"^) for QE{N,t), and hence so 
is the asymptotic behaviour of the small x splitting function, i.e. Pgg{x,t) ~ x~^''^^. So we see 
that the introduction of the running coupling has a dramatic effect on the singularity structure 
of the LO BFKL anomalous dimension, turning the cut into a series of poles, and changing the 
position of the rightmost singularity by a factor of ~ 0.4. This result of the pole in the anomalous 
dimension was previously proved in detail in [15] using numerical techniques and in the context 
of the coUinearly resummed NLO kernel, and also indicated here using an approximate analytical 
solution first suggested in [24]. However, in this paper I particularly stress the quantitative result 
of the huge modification of the naive LO BFKL anomalous dimension due to the running coupling 
contributions alone. This is apparent over a wide range of N, and in fig. 5. a I show the anomalous 
dimension as a function of A'^ for all values right of the leading singularity. As one sees, it is much 
closer to the simple as(Q'^)/N expression than to the naive BFKL result. 

Before going into more precise detail and more general situations there are two important 
points I should address. These are the choice of the scale of the running coupling in (2.11) as fc^ 
and the fact that the expansion of Xo(7) iii powers of 7 is not convergent over the whole range of 
the contour of integration. The former of these is the simpler, so firstly I shall address the choice 
of scale. It was known in [26] that the correct scale seemed as if it were really the symmetric 
choice {k — q)"^, but that k"^ could be used instead, leading to contribution to the NLO kernel 
which is proportional to (3q. In practice it is much easier to obtain analytic results using k"^ , and 
this /?o-dependent NLO term leads to a contribution to the Mellin transformation of the NLO 
kernel, Xi{l)i of the form \^o{xo{l) + Xo(7))- Including this in the integrand for the expression 
for Qe {N, t) at NLO (to be discussed in detail in section 4) leads to a multiplicative contribution 
of the form exp(|(ln(xo(7)) + A^o(7))) = f^°{l)- This can be expanded as a power series which at 
low orders is 



//3o(^) = 1 + 1.6O73 + 1.247^ - 0.1637^ + 1.157^^ + • ■ ■ . 
Including this additional factor in (2.23) modifies (2.32) to 



(2.34) 



V (PoA)t'^ 




(2.35) 
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For a given power of as{Q'^) these new contributions produce terms a power of ^o-^ up on the other 
terms and hence, not surprisingly, result in additional running coupling corrections to the gluon 
and anomalous dimension. However, the new terms in the series in powers of 7 do not start until 
third order and have rather small coefficients. The resulting change in the anomalous dimensions, 
both for general values of N and for the position of the leading pole is very minor. Therefore, the 
correction for my original "incorrect" choice of scale is very small. However, in principle it seems as 
though the factor just considered should really be taken as part of the LO result since it just gives 
running coupling corrections. I will adopt this convention and the LO anomalous dimensions and 
splitting functions presented in this paper will explicitly contain the corrections from this factor, 
and in fact the results already presented in fig. 3, fig. 4 and fig. 5 include these (very small) effects. 
In principle one could sum the corrections needed due to the simple choice of in the coupling, 
rather than {k — q)^, by including contributions induced in the kernel at NNLO and beyond. In 
practice, beyond NLO the change seems too tiny for one to be concerned. 

I should also comment on the limit of applicability of the analytic expression (2.30). As noted, 
it is obtained via a series expansion which is not valid over the whole contour of integration. This 
is reflected in the error of 0{A? /Q"^) we discovered for this expression but also in the fact that 
the overall magnitude of the Ajj(— l/(/3o-/V)) actually increases like n\ in general. This latter point 
means that the series in (2.30) is actually asymptotic. It turns out that it contains both infrared 
and ultraviolet renormalon contributions, and hence it must be truncated to obtain sensible results. 
The greatest accuracy may be obtained from (2.30) by truncating the series at order uq ~ t, the 
precise value depending on the size of the coefficients in the scries expansion. For the LO gluon 
these are small and one could use uq ~ 10, but from experience with other variables (sec later) 
and the desire to go down to IGcV^, i.e. t 4 — 5, in practice I always use uq = 5. (For 

the LO gluon the contribution from n = 6 ^ 10 is practically negligible.) Using the truncated 
expression for QE{N,t) in the manner already discussed, then results in an infinite series in agiQ"^) 
for ■^gg{N,t) which is convergent for any N right of the leading pole, but different to the real, 
divergent series beyond 6th order in as{Q'^)- 

It is vital to note that although the formal expression for the gluon, and hence anomalous 
dimension, as a power series in as(Q^) (2.30) contains infrared renormalons^, and hence has an 
ambiguity of ©(A^/Q^), the integral in (2.17), which properly defines the leading twist gluon and 
anomalous dimension, does exist and produces well-defined results. The ambiguity of 0{h? /Q"^) in 
(2.30) cancels with an ambiguity in the 0{h? /Q"^) correction to this power-series expansion which 

^ In unphysical regularization schemes, such as MS, the anomalous dimensions are not expected to 
contain renormalons (see section 3.4 of [27] for a discussion), these being confined to the coefficient functions 
relating the parton distributions to physical quantities. However, by regularizing via a finite Qo, and 
defining the gluon density as the bare density convoluted with the gluon Green's function we have implicitly 
chosen a more physically motivated factorization scheme which allows the presence of renormalons. 
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we discovered in the derivation of (2.30). The accuracy of the (truncated) analytic expression can 
be found by comparing with results obtained from evaluating (2.17) using numerical integration 
along the contour shown in fig. 1. For the gluon structure function for N to the right of the 
leading pole the analytic approximation to the anomalous dimension is found to be a fraction of a 
percent for t = 6, and falls like exp(—t). Strictly speaking there is an exp{—t) contribution from 
the correction to (2.30) (with the renormalon ambiguity removed) plus a correction due to 
the truncation. However, l/t"^ is similar to exp(— t) in the range considered. Hence, we have a 
power-like correction to the power series in as{Q^) obtained from the truncated expression which 
is completely well-defined. This illustrates that the presence of infrared renormalons in a physical 
quantity is not necessarily due to an inherent ambiguity in the quantity itself (due, for example, 
to the Landau pole in the coupling) as is commonly thought, but rather due to the impossibility 
of completely expressing the physical quantity as a power series in as{Q^) [28]. In truncating 
the power-series expansion in (2.30) I simply choose to split the expression for the gluon as some 
general function of N and into a perturbatively calculable part as a power-series in as{Q^) and 
a remainder which is approximately of order 0{K^ /Q"^). The point of truncation is then chosen 
empirically so as to make this remainder term as small as possible. This seems to be the way to 
obtain the most accurate analytic results. It is important to note that the remainder term, although 
power-suppressed, is not in any way higher twist, since it is obtained from the leading twist part 
of the solution to the BFKL equation. 

Having got these two points out of the way we can now begin to discuss the quantitative 
results of the running coupling BFKL equation. In order to investigate the real effect of the BFKL 
anomalous dimension on structure function evolution it is necessary to calculate the BFKL splitting 
function as a function of x. This is where an analytic expression for the anomalous dimension is 
particularly useful. A series of numerically obtained values of jgg{N,t) allows an approximate 
determination of P{x,t), but it is extremely difficult to be accurate, especially for the wildly 
oscillating functions of 1/N which do in fact make up GE{N,t). However, I now have an explicit 
series for jgg{N,t) in powers of as{Q^), obtained from the truncated expression for GE{N,t). The 
A''-dependent functions at each power of as{Q'^) become larger at small N as the series progresses, 
of course, and to reach small enough x more and more terms are needed. However, at a fixed value 
of N there is no such growth, and the same is therefore true for fixed x. Hence, one only needs to 
work to a finite order. Limiting oneself to x > 10~^ and t > 4.5, i.e. ^ IGeV^, the suppression 
of the A„(— 1/(^0-^)) is quite significant and seventh order in as{Q^) is easily sufficient. This 
results in a power-series contribution to the splitting function 

xPgf{^,a,{Q')) = a,{Q^) + at{Q^)(2A^ - 12.01^of + 9.206^o'^ - 9.60^3 



+ aliQ^) (2.08^ - 26.95^0^ + 13A.6P^^^ - 320.7 + 359.8^0^? - 148.8^o') 

+ aJ(Q2) f — ^ - 19.23|j + 78.94^0 1^ - 169.2^2 1[ + 199.8^3 |j - 122.9^9^ |- + m.72^U\ . 

(2.36) 
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This contribution to the sphtting function for f = 6 and is shown in fig. 6.a. Note that because 
of the truncation of QE{N,t), beyond 6th order the expression for PgJ^{^,as{Q'^)) is not what one 
would really get from the true power-series. In particular there are higher powers of ^ than strictly 
allowed. Nevertheless, it represents a very accurate approximation to the full result whereas the 
correct series would simply diverge. 

We also have to consider the power-suppressed contribution. Although this is only calculated 
numerically in A?^-space it is only a small correction of order 0.05% for jg^{N,t) a t t = 6, and 
can also be calculated for a wide variety of values of N and t without too much work. It can 
then be modelled by an analytic function which may easily be converted to x-space. Hence, 
I choose to calculate it for t = 4.5 (Q^ ^ IGeV^) and t = 6 (Q^ ^ 6GeV^) and N-values 
0.4,0.5,0.6,0.7,0.8,0.9,1.0,1.5,2,3,5,00. The lower t value is the lower limit at which we will 
trust this LO perturbative result, and for t above 6 the power-suppressed effect is very small. The 
N values go low enough to correspond safely to x > 0.00001 and are sufficient that very accurate 
modeling can be done. The values are fit to a function of the form 

ao exp(-6ot) + exp(-t) «n (^^^^j;^) " ^) ■ (2-37) 



n=l 



Introducing further degrees of freedom beyond this does not seem to change the results. This 
expression can then be trivially converted to x-space. Performing this procedure in the case of the 
power-suppressed contributions to the LO gluon anomalous dimension I obtain the explicit result 



(4.5)7 2! ■ V"s(4-5)/ 3! ' \a,{4.5) J 4! 
+ 3.00 \ ^-0.301' ^' ^ ^ 



(4.5)7 5! ■ \a^{4.5) J 6 

(2.38) 

This power-suppressed correction is shown along with the power-series part and the full LO splitting 
function in fig. 6. a. Although the power-suppressed contribution in x space turns out to be a larger 
fraction of the total than in AT-space, it still only makes a very small correction to the evolution. 
However, one notices than the logarithmic terms in (2.38) are such that it falls more quickly than 
{h? /Q"^), or alternatively, grows more quickly than this as falls. This may be due to the presence 
of a significant (A^/Q"^) term in practice. 

The full LO splitting function is shown in fig. 6.b along with the purely order as{Q^) contribu- 
tion and the naive BFKL splitting function. One sees that it is hugely suppressed compared with 
the naive LO BFKL splitting function, and is even lower than the 0{as{Q'^)) contribution for x 
between about 0.1 and 0.001. Finally I note that the LO running coupling BFKL equation has also 
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been calculated in [29], but numerically, with coupling scale equal to {k — qY, and with the coupling 
frozen below a particular scale and Qq taken to be a finite value. The results are displayed for high 
t (where my power-series is essentially exact) and despite the above differences seem to be in very 
good agreement with the results in [21] and this paper. The freezing of the coupling and the finite 
Qo introduce choice-dependent non-perturbative effects which become important at extremely low 
values of x, which in general become lower as Qo and the scale of freezing decrease. This seems 
to support the results obtained by my method of formally factorizing the non-perturbative effects 
into Qi{Qo,N) and extracting as much information as possible in an analytic model-independent 
manner. 

3. Higher Twist at Small x. 

In this section I will show that as far as the information from the BFKL equation is concerned 
calculable higher twist contributions are small. I will also suggest that some other powerlike cor- 
rections at small x may perhaps be less significant than often claimed. As a first point I note that 
it has been claimed that there are likely to be large infrared renormalon contributions to structure 
functions at small x [30]. As shown in the previous section for the case of the gluon both infrared 
and ultraviolet renormalons do show up in the solution to the BFKL equation if one insists upon 
trying to express results entirely in terms as a power-series in as{Q^) and uses the whole of (2.30) 
rather than truncating. Presumably these are an extension of the small x divergent contribution 
to the renormalons in [30]. However, these renormalons are circumvented if one considers the full 
solution to the Q^-dependent part of the BFKL equation. Precisely the same argument works for 
the case of real structure functions, as will be shown explicitly in section 5. This is not to say that 
there are not relatively large power-suppressed corrections to the (truncated) perturbative-series. 
We have already seen a non-negligible contribution to Pg^{x,Q^), and the power-suppressed con- 
tributions turn out to be larger for physical quantities. However, these contributions are calculable 
and unambiguous. Hence, solution of the BFKL equation, which provides results more general than 
a power series in as{Q^) avoids the renormalon ambiguity. This means that renormalons obtained 
from unresummed (in ln(l/a;)) calculations require not only a ln(l/a;) resummation but also the 
consideration of results beyond the power-series expansion. This implies they do not really tell us 
anything truely quantitative about power corrections in practice. 

Now let us consider genuine higher twist effects. Some of these are contained within the BFKL 
equation, since if Qq is allowed to be non-zero a series in powers of {Qq/Q^) is obtained which 
tells us about the higher twist contributions due to the intrinsic transverse momentum in the two- 
gluon operator. This is the only information, however, and we learn nothing about the other three 
contributions to next-to- leading twist (discussed, for example in [31]), in particular those due to the 
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four-gluon operator and hence possible saturation effects. However, it is possible to obtain some 
useful and interesting results. 

Let us first consider the fixed coupling BFKL equation. When solving (2.8) it is straightforward 
to also calculate the higher twist contributions by picking up the non-leading poles in 7. The easiest 
way to proceed is to obtain G{Q^ , x) by first taking the exact inverse Mellin transformation back 
to x-space by picking up the simple pole at iV = asXoil) resulting in 



1 r^+ioo 

xG{Q\x)^—j^ d7exp(7ln(QVQo))exp(e«sXo(7))- (3-1) 

This can now accurately be evaluated in the asymptotic small x limit using the saddle-point tech- 
nique, i.e. integrating along the contour determined by the condition (dxoin)/d'^) = which defines 
7o. At leading twist, > 3fie7 > 1, 70 = 1/^ £^nd Xo(7o) = 41n(2), leading to the usual power-like 
growth at small x. However, looking for the solutions to (dxo(7)/c?7) = for — 1 > 3fie7 > 0, i.e. 
examining the higher twist operator and its anomalous dimension, one finds 

7^^ = -0.425 ± 0.474? Xo(7^^) = -2.64 ± 2.393z. (3.2) 

Hence, the features of the saddle-point are completely different at next-to-leading twist. Not only 
are there complex conjugate saddle-points leading to an oscillatory behaviour, but the real part of 
Xo(7o^^) is negative rather than positive. Inserting (3.2) into (3.1) one obtains 

xG^^(Q^ x) oc x^-^^"^ cos(2.393a,0, (3-3) 

i.e. a valence-like gluon rather than one growing at small x. The corresponding higher twist 
splitting function has the same general behaviour as the gluon as x ^ 0. 

One can also find the splitting function by solving 1 — {q.s/N)xo{'^) as a power-series in {cXg/N) 
for the next-to-leading twist solution. This results in the explicit series 

'-\ /-\2 /-\3 /-\4 /-\5 /-\6 



^r(^./N) + l^{f -2if +2^ +4.4 = -29.2 = +80.2 



7 /-\8 /-\9 /_n10 



90.61 -298f^j +2084(^1 -6446 



/-xl2 ^_xl3 /_xl4 

+ 9157f^j + 20919f^j - 187924(^j +666008f^ 

/-\15 ^_n16 ^_n17 ^_\18 

-1.2xl0«(^i) +1.3X10«(^) +1.9X10'(^) -7.7X10' (9 
-1.7X10»(|)"-2.1X10'(|)"-2.0X10»(|)%.... 

(3.4) 
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which can be easily converted to x-space. The corresponding splitting function is plotted for ag = 
0.2 in fig. 7, and it clearly fits the expectation that xP^'^{x, ag) ~ x'^-^ cos(0.5^) as a; — > 0.^ Hence, 
although the first term in the series is the same as at leading twist, and implies a growth at small x, 
the summation of the series is extremely different, and the next-to-leading twist contributions from 
the BFKL equation are not only suppressed by {Qq/Q'^), but also become negligible at small x. This 
can also be shown to be true for the even higher twist contributions using the same techniques. 
This highlights the danger of using low order terms in the series for the splitting functions to 
estimate higher twist corrections, as in [31] - the summation of leading ln(l/x) terms may be very 
important, in this case of the two-gluon operator leads to a complete change of conclusion on the 
import of higher twist. Unfortunately, there is no knowledge at all of the corresponding series for 
the four-gluon operators. 

Given that the results from the fixed coupling BFKL equation were altered so dramatically 
at leading twist by the inclusion of the running coupling, we should see what happens at higher 
twist. As already mentioned, the higher twist contribution to the running coupling BFKL equation 
is given by 



X [\M-lHQl/^^)+Ml)/{PoN))d^gBiQlN), 



where the contour in the first integral has been moved to the left since the leading singularity at 
7 = is eliminated by the second integral. 

Let us consider first the case where t = ln((5^/A^) ^ = lii(Qo/^^)' which would be the case 
for deep inelastic scattering. Let us also, without justification for the moment, let the lower limit 
on the second integral be a constant, k ~ —1, so that we have factorization imposed. In this case 
we can evaluate the two integrals separately. Both the integrals can be calculated accurately using 
the saddle-point method. Thus, using the type of steps outlined in (4.1)-(4.5) of [20] one obtains 

exp(f\na,{q^)/N)dlnr 

(3.6) 

for the first integral and 

eicp(-f"^^^{as{q^)/N)d InqA 

\ Z.^ (3.7) 

[-x'oiin^siQD/Nw 



Unfortunately, because of large cancellations, the first 21 terms in the series for xPgg (x, as) are 
needed for x > 0.00001. 
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for the second. It can be verified numerically that these expressions are indeed good approximations 
to the precise results. Combining these we get the full next-to-leading twist gluon Green's function. 

^w(lS^7riMQ^)/N)d Inq^ 

1 r- (3-8) 

7ri^sm/m-x'oiiri^s{Q')/mn-x'oiiri^siQi)/Nw 

Hence, the anomalous dimension for the higher twist operator is simply that obtained for the fixed 
coupling, but with the coupling constant allowed to run with the scale, while the normalization 
is (roughly) the root of the fixed coupling normalization evaluated for as(Q^) multiplied by the 
same for Q!s(Qq). Hence, the result is much the same as for the fixed coupling case, with both the 
splitting function and the normalization decreasing and oscillating as a: ^ 0. 

It order to justify this conclusion it is only necessary to explain why we could assume the factor- 
ization. To do this we note that the saddle-point for the first integrand is at t = (l//3o-A^)xo(7o^"^(i)) 
and similarly for the second integrand with t ^ to. However, since t ^ Iq ^q'^ {to) is significantly 
to the right of ^^'^{t). The value of exp(— 7to + ^o(7)/(/3o-^)) along the real axis along with 
7 = 7o^^(to)) 7o^^(0 is shown in fig. 8. It is simple to rewrite (3.5) in the equivalent form 



g^'^' iQ^N) = — / -exp(7ln(gVA') -^o(7)/(^oiV))cf7 



/ 



exp(-7ln(Q2/A2) + Xo{j)/{(5oN))dj 



exp(-7ln(Q2/A2) +Xo(7)/(/3oiV)) 



(3.9) 

Using fig. 8, and remembering that the saddle-point integral for the first integral is parallel to the 
imaginary axis, and that the integrand very quickly decreases away from 7o^"^(i), we conclude that 
the value of the second integral in the second line of (3.9) is negligible compared with the first. 
Also noting from fig. 8 that there is little change if we alter the lower limit of the first integral in 
the second line to A; — 1, we obtain the factorization assumed above. Hence, in this t ^ to limit 
we find that we obtain factorization of the next-to-leading twist solution and that as for the fixed 
coupling case this is negligible as a; ^ 0. 

Even if to approaches t the results can be shown to be similar by numerical calculation. For 
example, in the extreme limit of t = to the first integral in the second line of (3.9) gives only 
half the saddle-point contribution, but one can check that the previously negligible second integral 
now gives a roughly equal contribution for all N. However, factorization is now clearly broken. 
Detailed numerical investigation shows that for to not much smaller than t we can write the higher 
twist contribution in the form {Qo/Q^)f{Q^,Qoi ^) where the total is a function of N which grows 
slowly with N, approaching a constant as A?^ — 0. This is consistent with the form cos(61n(l/a;)) 
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which we get for the factorized next-to-leading twist solution (the Mellin transformation of which 
is {N + a)/{{N + a)^ + 6^)), and certainly confirms that the gluon Green's function is falling as 

Therefore, the higher twist operators and their anomalous dimensions derived from either 
the fixed coupling or running coupling BFKL equation are negligible at small x, and for these 
higher twist contributions the use of the running coupling equation does not qualitatively change 
anything. However, we are currently not able to say anything about the contributions from the 
four-gluon operators, and hence about shadowing corrections etc., beyond relatively simple results, 
e.g. anomalous dimensions in the small x limit at LO in a^. There have been various suggestions 
that such shadowing corrections are large, but I feel that these estimates may well be severely 
exaggerated by the use of the approximation of this LO in anomalous dimension, and also by 
the fact that the even more restrictive double-leading logarithmic approximation is often used. 
This often seriously overestimates the size of the anomalous dimensions, coefHcient functions, and 
also the gluon distribution. I hope I have demonstrated that for the evolution of the higher twist 
two-gluon operator the LO-in-a^ double-leading-log approximations is indeed totally misleading. It 
is also interesting to note that a more complete calculation of the higher twist coefficient functions 
for the evolution of F2{x,Q^) due to the four-gluon operators [32] implies that the double leading 
log approximation is a vast overestimate. Even using very small values of the screening length 
{R = 2GeV~^ rather than the more usual R ~ lOGeV"^) and the very large LO GRV gluon 
distribution [33], it seems that the shadowing correction is almost negligible in the perturbative 
HERA range. Saturation effects will no doubt eventually set in for low enough x and Q^, but 
presently I feel the technology is not such as to predict where with any real accuracy. Certainly, 
resummations in ln(l/,T) tend to decrease the size of the gluon extracted from data, and this 
combined with the above considerations suggests a much smaller saturation effect, and total higher 
twist effect, than often supposed. Certainly the model-independent "rule of thumb" for strong 
saturation contributions that dF2{x,Q^)/dlYiQ'^ ~ Q'^a(x) and hence dln(F2{x,Q^))/dlnQ^ « 1 
is not even closely approached for any HERA data with > IGeV^. 

However, I note that in my examination of higher twist I have not examined the mixing between 
leading twist and higher twist operators or included any nonperturbative contributions due to, for 
example, the behaviour of the coupling constant at low scales. These two effects are related to 
each other. Such questions have been considered for toy models in [15] and [29], and numerically 
for the full LO running coupling BFKL equation [29]. These papers have considered the full 
anomalous dimension defined by dln{Q{Q^ , N))/dt, and the way this is affected by the higher twist 
corrections, rather than just dln{Q^'^ (Q"^ , N)) / dt considered above. They demonstrate that there 
are potentially serious modifications to the leading twist anomalous dimension due to the higher 
twist corrections introducing sensitivity to the form of the normalization factor Gi{Qq,N) which 
depends on the regular ization of the coupling at low scales and on the Qq dependence. Depending 
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on the assumptions about the nonperturbative physics, these contributions can be be important at 
extremely small x, generally changing the precise form of the power-like behaviour, and for more 
severe imposition of nonperturbative effects, i.e. letting them set in at higher scales, introducing a 
completely different asymptotic behaviour. Unfortunately, within the framework of my paper the 
formal divergence of Qi{Qo,N) makes a similar study impossible and, as mentioned at the end of 
the previous section, I simply have to appeal to these alternative results, in particular the smallness 
of X at which the power-suppressed modifications set in, in order to support the reliability of my 
more formal calculations. However, I also note that the smallness of the higher twist operators and 
their anomalous dimensions calculated in this section suggest that whilst these contributions from 
non-perturbative sources only set in at low or very small x indeed it seems perfectly possible 
that they will give a comparable, or even larger contribution at low x and low than the genuine 
higher twist contributions. 



4. NLO Corrections. 

In section 2 I demonstrated that using as(/c^) in the BFKL equation, as in (2.11), has a 
profound effect on the form of the solution both for the normalization and for the anomalous 
dimension. However, given the first conclusions regarding NLO corrections in the essentially fixed 
coupling case, it is particularly necessary to check that the results presented are not severely 
modified by the inclusion of the NLO kernel, i.e. that the perturbative calculations are stable. The 
NLO kernel was presented in [11] and the way in which to solve at NLO with a running coupling 
was presented in [14]. Writing the NLO equation as 

fik^Ql) = h{k\Ql) + (^^) ^{Ko{q',k') - a,{k^)K,{q\k''))f{q% (4.1) 

and using just the one-loop expression for the coupling'' leads to a 2nd order differential equation 
in 7-space 

'^hl.N) ^ _ 1 d(xo(7)/(7,A^)) _ ^^ u 2) 

This can be solved in a very similar way to LO, i.e. it factorizes into the same form as (2.17) with 
(5^-dependent part given by 

1 /-i+ioo 

Qh{N,t) = — I -eMlt-X,{j,N)/{PoN))dj. (4.3) 



^ Using the full NLO expression for the running coupling would lead to a huge degree of complication, 
and this has never been attempted. Since, so long as A is chosen appropriately, the one- and two-loop 
couplings are very similar, I do not imagine any major errors in the results below. 
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However, Xi{'y,N) is rather more complicated than the previous ^0(7)- It can still be expressed 
in the form 

X^ ij,N)= r XNLO (7, N)dj, (4.4) 

but now XNLoili ^) can be written as a power series in beginning at zeroth order with Xo(7)- 
As seen in [14], though here ignoring resummations in N, the explicit form is 

x.Loh, iv) = x.(7) - iv^44 + ^ (- (44) ' - A (^44) V ■ ■ ■ • 

Xo(7) Xo V VXo(7)/ VXo(7)/ / 

where the currently unknown NNLO contribution to the kernel, X2(7)) would also appear at order 
A^^ in principle. 

As already discussed in section 2 there is a contribution to xi (7) from the /3o-dependent terms 

induced by an "incorrect" choice of the scale for the coupling - fe^ rather than (k — q)^. Taking 

this contribution to the term in (4.5) which is linear in N, and combining with the LO expression 

we find the previously discussed result of only a minor change in the anomalous dimension and 

splitting function extracted. Hence, the choice of as{k'^) is reliable, and is easily corrected for. 

In this section I consider the rest of the NLO correction to the kernel, which is much larger, and 

henceforth I denote Xi(7) as the NLO kernel with the /3o-dependent part |;9o(xo(7)+Xo(7)) already 

extracted, and include the multiplicative factor f^°{'y) in the integrand in (4.3). This still leaves a 

decision as to precisely what I take "the NLO calculation" to mean. There are various possibilities. 

I could work at the level of the NLO correction to the kernel, and hence the BFKL equation, and 

solve (4.1), producing the infinite series in (4.5). Alternatively, I could truncate XNLoilj^) in 

(4.5) after the second term. However, doing this still leaves the question of whether to use the 

whole of exp(l/^o Ji (Xi(7)/Xo(7))<^7) or just expand it out to first order in P^^. 
2 

There are particular problems associated with all choices. If one solves using the full NLO 
corrected kernel then there is an infinite series in powers of N to consider in (4.5), which turns out to 
be important in practice (see below). Also, the gluon Green's function and anomalous dimensions 
obtained from this solution contain many subleading terms beyond just LO and NLO in ln(l/a;) 
(and running coupling type corrections to these), as is essentially obvious from looking at (4.1) - 
iteration of / leads the last term producing NNLO then NNNLO and so on. Hence, this choice is 
disguarded. If one instead truncates (4.5) at order N, one still generates a subset of higher order 
terms beyond those one wishes, though it is possible to proceed in this case at least. One can see 
the explicit form of the solution by sTibstituting the truncated (4.5) into (4.3) and proceeding as in 
section 2. The contribution to ^1(7, A^) coming from the second term, —N{xi{'j)/xo{"y)), leads to 
an expression of the same form as in (2.21), i.e. 

CXD 

Xi(7, A^) = Xo(7) - QiVln(7) - iVco - iV ^ c„7", (4.6) 

n=l 
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where the c„ may be calculated easily by performing a power-series expansion of the known functions 
of 7, i.e. 



c„7" = 0.4247 + 0.8057^ + 0.52l7^ + 2.2907^ + 1.2877^ + 2.9807^ + 



(4.7) 



n=l 



Hence, the integrand for G\;{N^Q'^) becomes 



7 



-{l-ciN)/{0oN)-lfl3o, 



/°(7)exp 7t - JWW^^^^ + Y.ianl^''^' - Ncn^)) 



Performing precisely the same type of manipulations as in section 2 results in the expression 
gUN,t) = -sin(^^lii^|p)r(-(l - ctN)/{PoN))e^p 



(4.8) 



X 



n=l m=l '- 



1 + i„(l/(^oiV)) 1 + Cmil/Po) - 1 



[1-ciN) 



where 



oo ^ ^ CO \ 

^1 V/?o J 



(4.9) 
(4.10) 



and the A{1/{PqN)) include the contributions from /^"{-y), i.e. are of the form in (2.35). The 
factoring of the terms independent of t then results in the expression 



OO oo 



1 + An{l/{PoN)) 



n=l m=l 

X ( 1 + Cm{l/Po) ) - 1 



j.—n—m A 



-jl-ciN) 



(4.11) 

There are two sources of corrections beyond NLO in ln(l/x), other than running coupling correc- 
tions, in (4.11). Firstly, Cn{l/^o) can be expanded as a power series in l/(^o)- Only the first term 
in this series is genuinely a NLO correction to the LO result. Terms of higher order lead to con- 
tributions to the anomalous dimensions which are beyond NLO in ln(l/a;) without compensating 
factors of /?o which would enable them to be interpreted as running coupling corrections. Secondly, 
when one expands terms of the form ((1 — ciN)/{PoN))^ which appear in the A„ in (4.11), one 
obtains a power-series of the form. 



(1 - ciN) 



1 - nciN 



nin — 1) , ,,,9 



(4.12) 



The second term in this series gives the NLO in ln(l/x) correction while the remainder give higher 
corrections without compensating powers of /3o- Therefore, both these power-series expansions, i.e 
of the Cn in powers of l/(^o)) and the A„ in powers of N should be stopped at first order in 
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or N, and the cross-terms coming from first-order in both expansions, which are of overall second 
order, should be eliminated to obtain truly NLO results.^ 

Ultimately I define NLO by appealing to the perturbative form of the gluon Green's func- 
tion and anomalous dimension produced and hence by choosing the NLO definition such that the 
Green's function does receive only corrections which are no more than one power of as{Q^) (with- 
out compensating factors of /?o) down on the leading order one. This means using an expression 
for the gluon Green's function of the form 



QUN, t) = t(i-<^)/(^oiv) (l + £ f; [(l + ^(l/(/3oiV))) (l + cJPo) - 



(4.13) 



where the Cn/Po are obtained by expanding the exponential expression exp(l//3o(/7 (xi(7)/xo(7) + 

2 

Q + co))dj), out to just first order in 1/Po. Implicitly there is also a factor of 
— sin ^ ) ^ ~ ('^^^^^^ ) e^P(~7B/(A-^) + cq/Po) which contributes to the normalization 
in (4.13). 

Now that we have this NLO expression for the gluon Green's function it is necessary to make 
one more decision regarding the definition of the anomalous dimension. This is obtained from 
^LO+NLO^jy^^^ = {d\n{gl;{N,t))/dt). However, strictly speaking, in order to obtain only NLO 
contributions to the anomalous dimension {Ql,{N,t))~^ in this expression should be expanded only 
to NLO. This leads to a formal problem already pointed out in section 6 of [20]. Using the whole 
of {Ql,{N,t))~^ in the expression for the anomalous dimension we notice that the position of the 
first zero is changed from that at LO, leading to a shift, in fact a decrease, in the leading pole for 
the anomalous dimension, and hence in the power of leading behaviour of the splitting function 
as a; — > 0. So the x ^ behaviour of the splitting function becomes Pgg{x) = exp(Ao^ — AA^). 
However, since AA is due to NLO corrections, the strict NLO expansion is just Pgg{x) = exp(AoO ~ 
AA^ exp(AoO- This definition does not explicitly retain the shift in the power-like behaviour, and 
also leads to the NLO correction ultimately becoming larger than the LO result. Hence, I choose to 
retain the whole of (^^(AT, t))~^ in the definition of the NLO anomalous dimension, thus obtaining 



Ignoring this requirement and using the whole of (4.11), it turns out that the resultant expression 
is very badly behaved - blowing up at large N. This is almost entirely due to the higher-order terms in 
the expansion of the A„. Using the full C„(l//3o) does not change things much in practice. This large N 
instability translates into huge corrections in the splitting function at large x. Presumably this instability 
at large N and x is cured if one resums the whole series in (4.5). Including just the 0{N^) term does seem 
to improve matters. 
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the full Pgg{x) = exp{Xo^ — AX$,) as x ^ 0, even though in practice the choice makes little difference 
at the values of x relevant to HERA. 

So now I can use (4.13) to determine analytic expressions for the NLO gluon Green's function 
and anomalous dimension. However, the formal definition again results in a divergent power series, 
and as at LO I really truncate the series in (4.13) at no = 5. This leaves the problem of calculating 
the power-suppressed corrections. In order to do this it is necessary to have an exact definition 
for QE{N,t) in the form of an inverse Mellin transformation, as in (4.3). This requires finding 
the integral expression which would lead to (4.13) if a power-series expansion of the integrand is 
performed. Unfortunately this is not that simple. The problem comes with the manner of treating 
the —ciN term in (4.6). In order to have the leading t(i-'=i^)/(^oJV) factor in (4.13), and hence 
obtain the correct expression for the 0{as{Q'^)) part of the anomalous dimension, it is necessary 
to keep —ciN\n{^) in the exponential in the integrand, giving a factor ^-^i/Po. Expanding out 
exp(— c; ln(7)/^o) to first order would lead to ln(t) contributions to the anomalous dimension. 
However, keeping the full ^"^'/^o factor results in the argument of the A„ being —(1 — ciN) / {^qN) 
as in (4.11). Hence, there is no simple way to generate only NLO corrections from this term. In 
order to obtain an expression equivalent to (4.13) I choose to effectively put the known factor of 
f(-ciN/{/3oN)) |-,y hand and to generate the derivatives of the A„ within the integral with respect 
to 7. 

In order to sec how to do this I consider the LO expressions (2.17) and (2.26). It is quite 
simple to generate the first part of (4.13). All one needs do is insert the series expansion 1 -|- 
l/(/^o) X^^i Cnl^ expanded to first order in 1/Po into the integral representation, i.e. 

G'/{N,t) = [ j-'/^^^^^-'f''ij)e^pU - 7^ EanT^-^^) fl + ^ (l/^o)c^7"^ V^' 

(4.14) 

where the integral is over the full, unspecified contour, and generates the t-independent factor 
sin ^ (fi^u) ^ ^ (PoN) ^ ' ^ ^ t-dependent parts explicitly in (4.13). On top of this one 
must also insert the t-'=iJV/(/3oAr) factor by hand. If one is also concerned with the A?^-dependent 
normalization it is probably most consistent to also multiply by the factor 

sin(7r(l - QiV)/(/3oiV))r(-(l - qAr)/(;3oiV)) exp(-(7E - c^N)l{fi^N)) 



sin(7r/(/3oiV))r(-l/(/3oiV)) 
in order to obtain the overall factor of 



(4.15) 



-sin((7r(l -QiV)/(;3oiV))r(-(l -qiV)/(/3oiV))exp(^^:|±^). (4.16) 

Generating the second part of (4.13) is rather more complicated. One has to somehow modify 
the integral representation so that the derivatives of the A„(— l/(^o-^)) are obtained. To see how 
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to do this we let 1/{(3qN) = z, in which case the equivalence of (2.23) and (2.30) (ignoring the 
divergence of the series) is 

f 7-^-1 exp Lt-zf^ anl^'^+A dj = - sm{Trz)r {-z)t' (^ + fl ^n(^)i""A(-z + n)) , (4.17) 

•^C' V n=l ' ^ n=3 ' 

where I have removed the trivial factor of exp(— 7e/(;5oA^)) from each side. Differentiating both 
sides with respect to z we obtain 

- ^ ln(7)7— 1 exp [^ft - z f;a„7'"+') d^i 

^ 7"""' E exp ( 7i - z a„7'"+' ) 

m=l ^ n=l ^ 

= ^-(-2) sm{Trz)T{-z)t^ + £ ^„(2)t""A(-2 + n)"] 

^ n=3 ^ 

- TT cot(7rz) sin(7rz)r(-z)t^ ( ^ + XI -4„(z)t""A(-z + n) J 

^ n=3 ^ 

- ln(t) sin(7rz)r(-z)t^ + A„(z)t-'^A(-z + ra)") 

^ n=3 ''^ 

- sin(;r.)r(-.)i^ (1 + f; A„(.)r-^^^^^^) 

- sin(7rz)r(-z)i^ (^f; ^^£lt-"A(-^ + n)^ . 

(4.18) 

The last term on each side are equivalent, and rearranging the rest we obtain an expression for a 
series containing the derivatives of the A„(z) - 



sin(7rz)r(-z)t^ (^f; yl„(z)t-" ^^Lj±^^ = ^ ln(7)7— ^ exp {^t - z , 



(00 
1 + ^ A„(z)t""A(-z + 
n=3 



(4.19) 



which using (4.17) becomes 
sin(7rz)r(-z)i- ff;A„(^)t-" ^^^7 + ^^ 

^n=3 



ln(7) - (*(-2;) - 7rcot(7rz) - Int) 



7-"-' exp (7* - ^ E d7. 



(4.20) 

Therefore, the right-hand-side of (4.20), multiplied by —ci/{(5o)t^~'^^^^^l^°^'>\ gives the second term 
in (4.13) with some t-independent normalization which should be multiplied by (4.16) to be con- 
sistent with the first term in the preceding paragraph. Thus, we have a prescription for the full 
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calculation at NLO which is equivalent to the series expansion in (4.13), i.e 



^0 [■ 



m=0 
oo 



d7, 



(4.21) 

and once again one should multiply by (4.16) to get the most suitable normalization. We can now 
insert the above expression into ^^^'^^^^{N,t) = {d\n{Q\,{N ,t)) / dt) and evaluate numerically in 
order to get the NLO anomalous dimension without recourse to the truncated series expansion. 

We are now in a position to solve for the anomalous dimension and splitting function at NLO. 
Unlike the case of fixed coupling, or the simplistic results of the saddle-point evaluation, the NLO 
corrections to the LO anomalous dimension are under control. This is simply illustrated by the 
positions of the leading pole in the anomalous dimensions which are shown in fig. 4, and one can see 
that they change from about 0.25 for jgg(N, t) at LO to 0.17 at NLO, and that the Q^-dcpendence 
reduces a little. However, as already noted at LO, the value of the intercepts has little to do with 
physics at HERA - the power-like behaviour only really settling down for lower x, and this is even 
more true at NLO. Being more particular one notices that the anomalous dimension 'jgg{N, t) over 
a wide range of shows only a relatively small change going from LO to NLO. This is shown in 
fig. 5.b where the part of the NLO anomalous dimension at first order in as(Q^), i.e —0.935as{Q^), 
is not included, since this should properly be included at LO in a combined leading order in as(Q^) 
and as(Q^) ln(l/.T) expansion scheme. Alternative definitions of NLO lead to very similar results 
except at very high values of A^, where less sophisticated definitions lead to blowing up at large A^, 
as already mentioned. For this case of the gluon structure function the NLO correction is negative 
except for very large A^. I should also note that the power-like correction to the purely analytic 
result is a larger proportion of the NLO correction than of the LO contribution, but would still be 
almost impossible to spot if shown on fig. 5.b. The correction to the analytic value for the intercept 
is about 7% at t = 6 however. 

One can also make the transformation to x-space and calculate the NLO corrected splitting 
function. Unfortunately, due to the increase in size of the c„ coefficients compared to the a„ 
(particularly the absence of zeros) and also to the factors of n invoked by differentiating the A„ 
in (4.11) the power-scries in ctgiQ"^) is much less convergent than at LO. In order to obtain an 
expression which is reliable down to x = 0.00001 at = IGeV^ it is necessary to go to 20th order 
in as{Q^)- Hence we can write the NLO correction to the splitting function as 

19 r rirna x / j-jji 7jn — m — \ \ 

xPj^^O(e,a,(Q2)) = a^(g2^^ ^ ^n^Q2^UJ_J2^^^ + K^,-pn^(l-x)\, (4.22) 

n=l m=0 ^ ■ ^ 
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where because we have truncated the series for the gluon structure function rrimax can be greater 
than the naive expectation of rrimax = n — 1. The coefficients for the series are shown in table 1. 
If one is only concerned with x > 0.0001 or > 4GeV^ then the series can be truncated at about 
12th order. 

As at LO we also have to model the N dependence of the power-suppressed correction by an 
analytic function. Fortunately, exactly the same type of function is sufficient and we obtain the 
power-suppressed NLO correction to the splitting function of the form 

88 1 21 

-2.86exp(-1.02Wl-x) + expH)(l3.59(£=§^)' _29.6l(^)- { 
+ 39.76 

- 4.479 

The full NLO correction xP^g^'^{x) and its power-series and power-suppressed contributions 
are shown in fig. 9. a, where the relatively unimportant terms oc (5(1 — x) are absent. As at LO the 
power-suppressed correction is proportionally much larger in x-space than in moment space and 
certainly needs to be considered at t = 6 and below. Also as at LO it tends to oppose the form of 
the power-series expression, hence reducing the total NLO correction. The powers of in (4.23) 
are slightly smaller than for LO, and hence the power-suppressed correction does not fall quite so 
quickly with . 

The total NLO splitting function, i.e LO plus the NLO correction, is shown for t = 6 in 
fig. 9.b, where the contributions oc 8{l — x) both from the 0{as{Q'^)) part and the running coupling 
corrections to this are absent. The latter of these is a very small contribution. The NLO corrected 
splitting function is clearly not qualitatively different from that at LO, though it is quite a lot 
smaller at small x. Hence it seems as though by including the infinite series of running coupling 
corrections the perturbative expansion of the BFKL splitting function has been stabilized. However, 
the real importance of the NLO corrections as far as physics is concerned is the effect they have on 
the evolution of the gluon structure function. This is demonstrated in fig. 10 where the evolution 
of a suitable model for the structure function G(a;, Q^), i.e. (1 — is shown both for the LO 

running coupling splitting function, and for the NLO corrected one (all 5{1 — x) contributions other 
that at first order in as{Q^) one are included). Also shown is the evolution due just to the double- 
leading-log term P{x) = as(Q'^)/x. As one sees, at this (fairly low) value of t, i.e. ~ 6GeV^, 
the evolution driven by the LO splitting function is very similar to that from the double-leading-log 
contribution, and is even slightly smaller for x from 0.007 to 0.00001, corresponding to the dip in 
the splitting function seen in fig. 6. Below this the growth of the splitting function increases the 
evolution above the double- leading- log result. One also sees that the effect of the NLO corrections is 



a«(4.5) 
a«(4.5) 



2! 



- 33.765 



2.16 



««(4.5)y 3! 

2.63 ^6 



+ 16.89 



as (4.5) 



(^) 



4! 



f4.23) 
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certainly significant, and increases relatively with falling x, but it is clearly a correction rather than 
the complete change in qualitative behaviour induced by the NLO corrections without resummation 
of running coupling effects. 

A further way often used to investigate the perturbative stability of a fixed order perturbative 
calculation is to investigate the renormalization-scale dependence. This is often used fallaciously, 
e.g. if one calculates Pgg {ag , x) to NLO in the standard perturbative expansion and then investi- 
gates variation of renormalization scales one will never notice the influence of the terms at higher 
orders in which are also of higher order in ln(l/x). This is symptomatic of the fact that the 
expansion purely in powers of ag is not really a correct expansion scheme for splitting functions 
(for a full discussion see [8]). However, once we have performed a resummation of large logarithms, 
as here, renormalization-scale variation should be more reliable. The renormalization scheme de- 
pendence may be investigated by letting 

a,(g2) ^ a^ikQ^) + l3oln{k)alikQ^) (4.24) 

and in the LO part of the splitting function expanding out to first order in ln(A;), whilst in the NLO 
part using only the zeroth order, i.e. just letting as(Q^) — > agikQ"^). In this case we must also use 
a similar procedure for the power-suppressed corrections, i.e. these are really of the form (A^//x|j) 
rather than {JS? /Q"^). The results for k = 0.5 and k = 2 are shown in fig. 11 for ^ 6GeV^ As 
with the NLO corrections to LO the variation is significant but leads only to a correction rather than 
a qualitative change. This implies that the series expansion is stable, if not as rapidly converging 
as one might ideally hope for. 

Hence, the NLO corrections to the running coupling BFKL derived splitting function are well 
under control, both in terms of the asymptotic power-like behaviour of the splitting functions 
and in terms of the evolution in the range currently accessible to experiments. For deep-inelastic 
scattering, or indeed any process where there is factorization of the infrared physics into the input 
parton distributions, e.g. Drell-Yan scattering in proton-proton collisions, no further resummation 
is necessary, or even useful, beyond the running coupling corrections. This in distinct contrast to 
the case where both ends of the gluon ladder arc associated with a hard scale. In this case the 
conventional BFKL expansion is fundamentally flawed due to progressively higher order poles at 
7 = and 7=1 (corresponding to large logs in the ratios of the two scales k\ and as shown in 
[13]. These large order poles need to be resummed, and without this resummation calculations are 
badly behaved over the whole range of A'" (in fact explicit calculation shows that this is particularly 
the case at large A'^). In the case of deep inelastic scattering the coUinear factorization procedure 
automatically orders the poles at 7 = correctly, and the above problem shows up in high order 
poles at 7 = 1 only. The anomalous dimension is totally dominated by the region very close to 7 = 0, 
as this paper shows, and is very insensitive to effects at 7 = 1. Including the type of resummation in 
[13] [14] alters results from the NLO corrected case by only a very small amount, and is likely to be 
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no more influential than the remaining NNLO effects for which it does not account. Resummation 
of poles near 7 = 1 would be essential if one attempted to obtain information about the input form 
of the gluon, i.e. Qi{Qq,N). However, as well as the fact that Qq is an essentially nonperturbative 
scale, this type of calculation, along with the whole subject of single-scale processes, is also plagued 
by the infrared ambiguity problem caused by behaviour of the coupling at low scales. A discussion 
of such issues can be found in [15] and [29]. 

I close this section by noting that although the above results all look promising it is important 
to realize that they are all in a sense ambiguous because they deal with a particular way of defining 
the gluon parton distribution, which is a factorization scheme-dependent quantity. In this paper it 
is defined in a manner which is natural from the point of view of the solution of the BFKL equation, 
and which one may think of as perhaps a good "physical" definition of the gluon. However, it is 
very difi'crcnt from, for example, the gluon defined in the MS scheme. In order to investigate the 
real success of the approach in this paper it is necessary to look at the results for the real physical 
quantities - the structure functions. 



5. Small X Structure Functions. 

One may define a real structure function by a simple extension of the above methods, i.e. by 
including a hard scattering cross section at the top of the gluon ladder. This modifies (2.4) to 

J^i{Q\N)=a, —aiJk'/Q^)f{N,k^Ql)gB{N,Ql), (5.1) 

where o'i^gik'^ /Q'^) is the cross-section for scattering of a virtual photon from a gluon with transverse 
momentum k'^ . For the case of the longitudinal structure function this cross-section is well defined 
even in the limit fc^ 0, but for J^2i^:Q^) the cross-section diverges like ln{Q^/k^) as k^ ^ 
(for details sec [34]). This demonstrates that for Tl{x,Q'^) the solution in the leading limit 
factorizcs neatly into the gluon distribution and a multiplicative coefficient function, while for 
J^2iN,Q'^) there is interference at this order between the coefficient function and the result of 
solving the evolution equation including the anomalous dimension 0^7^^ (as, N). In this latter case 
it is simplest instead to differentiate with respect to ln{Q^) obtaining 

where '^'^''^iI^q2^ is finite as k^ 0. In this case, if we work in an DIS-type scheme, i.e. one in 
which the quark-gluon coefficient function vanishes beyond zeroth order, there is a simple factor- 
ization between the anomalous dimension asjqg{as, N) and the gluon distribution.^ 



^ Note that in this article I ignore the mixing with the quark input distribution in general for simplicity. 
However, it does implicitly appear in the NLO correction to the kernel, i.e. it is the NLO correction to the 
anomalous dimension eigenvalue rather than to 7gg which I use since this is the quantity directly calculated 
in [11] [12]. The contribution to this due to the quark mixing is very small in practice. 
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In order to progress it is first necessary to consider the overall factor of in the above 
expressions, and particularly its scale. One might think that it should be as{k^), and thus appear 
within the integrals with respect to k^. However, this could only come about due to double 
counting of diagrams, since the resummation of bubble diagrams required to make this equal to 
as{k'^) has already been performed in defining the coupling in the BFKL equation as as{k^). 
is the only remaining scale, so it must be the scale of this coupling. One can also justify this 
by considering the fact that there is a NLO correction to the input of the BFKL equation of the 
form — /Jqcks ln((5o/MR)'^(^^ ~ Qo) (coming from bubbles in a gluon propagator). Introducing this 
into calculations leads to multiplying each result by a factor (1 — /Jqcks lii(Qo/A*fl))- This splits into 
—Potts ln{Q^ / fi"^) + PoCXs ln(Q^/<5o)5 ^^'^ latter term is an infrared divergence which contributes 
to the one-loop gluon-gluon splitting function while the second goes into making the overall factor 
of as have renormalization scale Q^. 

Now removing the overall factor of as{Q^) (or in fact the normalization factor as{Q'^)Nf j (Stt)) 
from (5.1), and taking the Mellin transformation with respect to {Q'^ / h?) leads to the simple 
expression 

J^,(7,iV) = /ii,,(7)^(7,iV). (5.3) 

Thus we may solve for JTj [N^ t) in exactly the same way as for Q {N, t) , obtaining exactly the same 
divergent Q^-independent part and a Q^-dependent part given by solving 

:FEAN,t) = — I ^^/''(7)exp(7t-Xo(7)/(^oiV))d7- (5-4) 

This may be evaluated numerically, using the same contour as for the gluon, or in order to find the 
power-series solution we may proceed as with the gluon structure function by expanding the hi^g{^) 
(which were calculated in [34]) as a power series about 7 = 0. For the two cases we discussed above 
we have 

/iL,s(7)/^°(7) = 1 - 0-337 + 2.137^ + 0.677^ + 2.587^ + 2.997^ + 1.927^ + ■■■, (5.5) 

and 

^2,9(7)/^° (7) = 1 + 2.177 + 2.307^ + 6.677^ + 7.057^ + 12.927^ + 15.477^ + • ■ ■ . (5.6) 

It seems natural to absorb the (in some sense) NLO corrections from (7) into the contributions 
from the hi^g^-j) since they are of exactly the same form, whereas the other NLO corrections have 
inverse powers of /3o- Following the same steps as in section 2.2 then results in an expression 

:FeAN, t) = i^/^'^"^) (1 + B,,Ay{PoN))t-^AA-l/{PoN))) , (5.7) 

^ n=l ' 
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where the Sj^„(l/(^o-^)) are now determined not only by the power series in 7 obtained from the 
expansion of ^0(7), but also from the expansion of hi^g(j). In particular they now contain parts 
at zeroth order in l/{PoN). 

Using these results it is now a simple matter to derive the longitudinal gluon coefficient function 
at leading powers of ln(l/a;) plus running coupling corrections and similarly for the quark-gluon 
anomalous dimension, i.e. 



CL,g{a,{Q^),N) 



2, ^^^_ a,{Q^)NfJ^E,L{N,t) 
Stt gE{N,t) ' 



(5.8) 



with obvious generalization to ^qg{cis{Q'^)-,N). These moment space expressions may easily be 
converted to x-spacc. Truncating the series for the structure functions and the gluon at no = 5 
results in the perturbative series for xCL,g{ois{Q^),x), 



xCL,g{as{Q^),x) = 



ocs{Q^)Nf 



Stt 



(5(1 -x)- 0.33a, (g2) + 2.13a2(g2) _ 



d c3 c2 



+ a^(Q2)(8.4li^ - 54.45^0^ + 125.20^^ - 121.2^o'C + 42.0^o') 



4! 
-0.89 



3! 



7.76|^ _ 27.53;3o^ + 49.48/52^ - 44.59/3^^ + 15.77^^^) 



4! 



3! 



2! 



3! 



2! 



+ aI(Q2)(_i_ _ 33.4l|- + 164.8/3o|j - 419.3;3o'|r + 577.2^o'|r - 404.9^o'|7 + 112.9^o'C) 
+ - 72.27|^ + 335.7/30 ^ 



7! 



6 

61 



838.2^9^ + 1210^^^ - 1004^^^ 



4! 



3! 



+ 441.7/?; 



e 

°2! 



79.05^0'?) 



(5.9) 

However, as for the gluon splitting function we have to calculate the power-suppressed correction 
by evaluating the inverse Mcllin transformations numerically. This is done in precisely the same 
way as for the gluon, and results in the correction to xCL^g{as{Q^)-,x) of the form 



37r 



(-1.168 - 0.482t + 0.1106) exp(-t)5(l - x) + exp(-t) 

0.875 



(-4.685(^ 



+ 34.25 
- 17.94 



Qs(t) 
a.(4.5) 

as{t) 
a, (4.5) 



i - 59.47 



1.37^4 

4[ 



+ 3.365 



Qs(t) 
as (4.5) 

as{t) 
as{4.b) 



0.074 ^2 

2[ 

1.77^5 
5[ ~ 



+ 45.81 



(4.5) 

Ois{t) 



-3.026 



as(4.5) 



0.2942 



U.(4-5) 



0.78 ^3 

3[ 

1.78 ^6 \n 

6" 



(5.10) 

where in this case it was necessary to model the N ^ 00, i.e. the S{1 — x) part with a slightly more 
complicated form than previously. Both expressions have been shown in a form which is sufficient 
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for Q2 > iGeV^ and x > 0.00001. The full xCL,g{x,t) is shown in fig. 12 a. along with the two 
contributions above. Note that the S{1 — x) term at 0{as{Q^)) in the power-series is obtained 
from the inverse Mellin transformation of the limit as A?" — > of the full 0{as{Q'^)) coefficient 
function and in the figure we replace it by the full 0{as{Q'^)) contribution, 6a;^(l — x), for ease of 
presentation (it not being easy to represent the normalization of the (5(1 — x) term). The 5(1 — x) 
term is simply missing from the power-suppressed part, though this is insignificant. We see that 
the power-suppressed contribution is now a much larger fraction of the total than for the gluon, 
though it docs not increase as quickly with falling Q^. In fig. 12 b. we show xCL^g{x,t) along 
with the 0{as{Q'^)) contribution and with the naive LO BFKL result in this factorization scheme, 
which grows far more quickly than the resummed result. 

Similarly we can calculate the perturbative series xPqg{as{Q^),x), 



xPqg{as{Q'^),x) 



37r 



6{1 -x) + 2.17a«(Q2) + 2.30^2 (g2) _ 

e ..... 



+ a3(Q2)(5.07l^ - 15.21^0^ + lO.M^o') + octiQ^){8.80^ - 47.50;3o|^ + 81.02^^^ - 42.30^3) 
+ af(Q2)(l8.88^ - 156.7;9of +478.0^o'C " 620.4^o'C + 280.3^9^) 



4! 

J6 



2! 



+ 



4 

4T 

e 



- 44.15|j + 159.9p'o^ - 293.4^2 ll- + 269.7^o'|l- ' ^^-^^^o^) 



+ al{Q^){~^ - 86.53^ + 385.6^o 5, 



i3 

3[ 



17.15 
~8! 



t2 

2[ 
3! 



234.6^ + 1354/3o|^ - 4263;9o'S + 7882.9^30'^ - 8519^o'^ 



t2 

764.0/3^^-^203. 

e 



6! 



5! 



4! 



3! 



+ 4962^^1^ -1199^o'0 



3.97 g^" _ 51.57 e 
1|"T0! "^9! 



+ 269.5|- - 647.5/3o|r + 258.8/?o'|7 + 2451/3o'|7 



e 

7! 



6! 



5! 



6962^^^ + 8473^o'§ - 5145^o'f + 1259^o'e) 



and we have a power-suppressed contribution to xPqg{as{Q^),x) of the form 
MQ')Nf\,,,,_, , ,,,,,,,, , a,{t) 



37r 



12.86 exp(-1.52U)(5(l - x) + exp(-t) -14.31 



2.695 



-41.14 



+ 1.85 



a.(4.5) 
a. (4.5) 



3.03 ^2 

2] 

3.695^5 
5[ 



+ 25.34 



0.1693 



as(4.5) 
a. (4.5) 



3.20 ^3 



3! 

3.80 ^6 



a. (4.5) 
9.096 



+ 36.297 

as{t) 



(5.11) 



2.93 



a«(4.5) 



3.44 ^4 



as (4.5) 
4! 



(5.12) 
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The full xPqg{as{Q'^), x) is shown in fig. 13 a. along with the two contributions above. As with 
xCL,g{x,t) the 5(1 — x) term at 0{as{Q^)) in the power-series is replaced by the full 0{as{Q'^)) 
contribution which is 1.5x{x^ + (1 — x)^). Again the 6{1 — x) term is missing from the power- 
suppressed part, and again this is insignificant. In this case the power-suppressed part is tiny at 
t = 6, though from the large powers of as{Q^) in (5.12) we see that it grows very quickly at lower 
Q^. In fig. 13 b. we show xPqg{x,t) along with the 0{as{Q'^)) contribution and with the naive LO 
BFKL result in this factorization scheme, which again grows far more quickly than the resummed 
result. 

These above results, along with the LO gluon splitting function, allow for a LO in ln(l/x) (with 
running coupling corrections) calculation and analysis of structure functions. In previous papers 
[8] I have strongly warned against the use of factorization-scheme dependent splitting functions 
and coefficient functions within the ln(l/x) expansion. It is still true that it is always possible to 
make huge redefinitions of the unphysical parton distributions by factorization scheme changes at 
a given order (or even at all orders) but the changes invoked by transfer between the commonly 
used schemes are diminished somewhat by the reduction of the size of the splitting functions and 
coefficient functions by the inclusion of the running coupling effects. It is also true that many 
of the changes invoked by factorization scheme changes are themselves due to running coupling 
effects, and the resummation of these stabilizes the whole procedure a great deal. Hence, it is now 
possible to work in terms of these unphysical quantities if one wishes, without potential disasters, 
as long as the ordering of the expressions is done with particular care. Nevertheless, it is still very 
convenient in some ways to eliminate the partons completely and work directly in terms of the 
structure functions Fl[x,Q^) and F2{x,Q'^) and the physical anomalous dimensions [22]. In fact 
we can easily argue a case for improved stability. At LO the longitudinal coefficient function is 
positive and quite large at small x, and hence Fl{x,Q'^) will be enhanced compared to the gluon 
at small x. At NLO the gluon evolution is smaller than at LO. Hence, evolving down from a 
given gluon at very high (where everything is simpler and more reliable) the NLO gluon will be 
larger at small than the LO gluon. However, we expect the NLO corrections to CL,g{x,Q^) to 
be negative, and thus counteract this increase in the NLO gluon in the calculation of Fl{x,Q'^). 
Hence Fl{x,Q'^) is (probably) a more stable perturbative quantity at small x than G{x,Q'^). 
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The physical anomalous dimension which is most closely related to the gluon anomalous di- 
mension is 



Ignoring the mixing with the quark sector this is given in terms of the parton related quantities by 



where I will use the convention of ignoring the overall power of as{Q^) in the coefficient function 
which would just result in a single contribution of —(3oal{Q^) to (5.14). Using the LO jgg{N,t) 
plus running coupling corrections, and similarly for CL,g{N,t) we see that the latter gives entirely 
running coupling corrections, and the total is the LO 'ygg{N,t) with an extended set of running 
coupling corrections. This total expression could be calculated from the 'ygg{N,t) and CL,g{N,t) 
already calculated, but part of the advantage in using physical anomalous dimensions is that it 
reduces the number of perturbative quantities governing the structure function evolution, i.e. the 
four splitting functions and four coefficient functions used to define F2{x,Q^) and Fl(x,Q^) are 
reduced to four truly independent physical splitting functions. Hence, we notice that using (5.4) 
for the longitudinal structure function we can calculate VLL{N,t) and PLL{x,t) directly, rather 
than from (5.14). Of course, the two definitions arc equivalent, but the latter allows a single 
power-suppressed correction to be calculated rather than having to combine those for jgg{N,t) 
and CL,g{N,t) and thus the potential error is minimized. The asymptotic power-like behaviour 
for P[^£{x,t) is not identical to that of Pg^{x,t) and is shown in fig. 4. The difference is only 
relatively minor, but one sees that the power-like growth for Fl{x,Q^) is slightly smaller than for 
the gluon, and is also slightly less Q^-dependent. The result for the LO in ln(l/x) power-series 
solution xP^£{as{Q^),x) is unfortunately a little less convergent than the previous LO quantities, 
due to large coefficients generated in taking the derivative with respect to t of the expression for 
TL{N,t) (or of CL,g{N,t)). Hence, in order to obtain an expression which is sufficiently accurate 
for Q2 > iGeV^ and x > 0.00001 we need to go to about 12th order. This results in the explicit 



rLL{N,t) 



dln{TL{N,t)) 
dt 



(5.13) 



VLL{N,t)=-igg{N,t) + 



d\n{CL,g{N,t)) 
dt 



(5.14) 
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expression 



11.29/30^ + 12 



+ a^(Q2)(2.4^ - 11.29/3o^ + 12.94^o'? - 4-02^3) 



+ a^(Q2) (0.121^0 1^ + 37.85^2 1[ - 99.88^o'C + 61.92^^) 

+ «'(Q')(2|^ - 75.14/3o|f + 454.7;3o'f " W3Ap^^ + lOll^o'^ - 358.8;3o') 
+ al{Q^){~^ - 13.94^ + 23.68^0 1^ - 39.48/32^ + 121.9^o'f - 155.2^^^ + 61.14^o'C) 
+ a^(Q2)(_i6.91^ + 348.8/3o^ - 2087^2 1^ + 5522^3fi - 7305;9o^^ + 4754^o'f - 1215^o'0 
+ - IIS-^I, + 1173^0 1^ - 5052/32|^ + 12044^33 1_ _ i7444;34|_ 

+ 15528^o'f - 7859^o'f + 1728^o'C) 

+ - + ''-'4 - ^''-'^"l + ''''^°|[ + - 24427^0^1 

+ 61280^o'^ - 72753^0' f + 42720;3o'f - 9998/3o'e) 

+ ""(Q')(-^IiT + ^"^^-Olo! ~ 2988^0 1^ + 5290^0' |r - 22253^0' |r - 135896/3^ |j 

+ 321404/3^1^ - 425485/3^^ + 330620/3;^ ||- - 141370/3^ + 25747/3^^) 

2 82 ^13 1417^12 til eio d9 t8 

+ 295058/3o^ - 538834/3^^ + 658339/3^-^ - 499685/3^^ + 211914/3^^ - 38311/3^^). 
The power-suppressed correction is calculated in the usual manner and is of the form 
36.57exp(-1750^(l-x) + exp(-0(4.626(J;2^)"''^37.84(J||)"■% 

Va.(4.5)y 2! Va.(4.5)y 3! Va.(4.5)y 4! 

- 3.136(^1 " "i! + 0.1706f ^) "^ "^l . 

The anomalous dimension r;^2(A/", t) is plotted in fig. 14.a. Until AT is very small it is similar 
to jg^{N,t) and both are close to the common as{Q'^)/N contribution, though r££(Ar, t) is a 
little larger at large AT. However, at lower AT r££(A'^, t) dips below the others before eventually 
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(5.15) 



rising above as{Q'^)/N but staying below jg^{N,t). Clearly the effect of tlie additional coefficient 
function, and hence additional running coupling corrections, is to make TLL{N,t) dip significantly 
below the 0{as{Q^)) contribution as{Q^)/N for a region and to reduce the value of the intercept 
compared to the gluon structure function. The effective splitting function xPl^£{x,t) is shown 
in fig. 15. In fig. 15. a we see that the power-suppressed contribution is larger for xPf'f (a;,t) 
than it was for xP^^{x,t). In fig. 15. b we see the outcome of the comparison of the anomalous 
dimensions for Fl and the gluon. xPl^£{x,t) starts a little higher at a; = and the dip below the 
0{as{Q^)) part is considerably more pronounced than for xPgg^{x,t). Also, going to x ^ 10~^, 
we see that the splitting function dips again, showing that the subleading poles in the anomalous 
dimension may have large residues compared to the leading pole, and that the increase in xPl^£{x) 
with decreasing x is not monotonic. This corresponds to the significant fall of rLL{N,t) below 
(^s{Q^)/N at A/" ~ 0.6. The eventual rise of TLL{N,t) guarantees that the splitting function will 
eventually rise again with the calculated intercept, i.e. like at even smaller x. However, 

for t = 6 this asymptotic power behaviour does not set in until x < 10~^*^ and in the region of 
X ^ 10"^ xPff (x) even becomes slightly negative. For higher t even smaller x is required, e.g. 
t = 8 (Q^ ^ 30GeV^) needs x to become as low as 10~^^ before the power-like behaviour sets in, 
though the size of the dip before this is smaller than for t = 6. This illustrates very clearly that as 
far as phenomenology at HERA, or any foreseeable collider, is concerned the value of the intercept 
for the anomalous dimension is simply not relevant to the evolution of structure functions. Indeed, 
it is very possible that before the power-like behaviour has set it unitarization effects have already 
become important. For collider phenomenology it is the splitting functions over the relevant x and 
range which one needs, and this requires the sort of detailed calculation in this paper. 
One can follow exactly the same procedure for the other important physical anomalous dimen- 
sion defined by 

^^^ = r.UQ',N)^UN,Q% (5.17) 

simply by using the LO expressions for '^''^^[^^i ^ and J^l (A", t) . The power-like behaviour as x — 
is governed by the poles in !F{N,t) as in the previous case, so the position of the intercepts is 
identical. The power-series expression requires the first 10 powers in order to be valid over the 
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required range of x and Q^, so I write it as 



xPif (a,(Q^),x) = 



(5(1 -x) + 2.5a,(Q2) + al{Q%(, - 0.167^o) 



+ «'(<3')(f - 12.72^0? + 12.0^2) + a^(g') (7.007^ - 41.41^of + 61.42^o'C - 26.82^3) 



3! 



2! 



+ a^(Q2)(5 _ 52.95/3o^ + 253.0;3o'f " ^^^-WU + 238.32;3o') 

+ - 87.30|j + 409.7^30 1^ - 773.3^o'|7 + 621.7^o'|7 - ^^^-^^U) 

+ - 117.8|j + 591.4^0 1^ - 1701^o'|r + 2792^3|r - 2315^o'|7 + 741.5^o'0 

—4 38(1 

+ a?(Q^)(^| + 52.4li^ + 70.68^o| 
+ 6307^o'f -1596^oC) 



1954/32 +6623/3^ 



4! 



9500/3^ 



3! 



10 



4.64^ 
^10! 



95 79 f ^ f ^ - 

^9! + ^^^-^ll " ^^^^-^^^^^ + 450.0/32^ + 9410^3:^ 



o|^ + 9410/33| 



- 26327^9^^ + 33805;3o5^ - 21743^0^^ + 5614^«^) 



4! 

7.478 



2! 

115.7^10 , ^^^^e e , „.„..2C' 



- 35038/3^ 



^ + 104536/3^-^ - 127198/3^ 



/3o 10! 
5! ' ""^^4! 
The power-suppressed correction is 



+ 765.7^ - 3293/3o^ + 8687/3o^^ - 5511(3^^ 



9! 



- + 74260^8 _ 



8! 7! 



6! 



17101/3^^) 



(5.18) 



(3.558 + 0.4216t - 0.1542^^) exp(-t)(5(l - x) + exp{-t) ( 72.17 

1.66 / \ 2.66 



78.03 



+ 13.50 



- 0.5141 



a. (4.5) 

a. (4.5) 

/ asjt) 
U.(4-5) 



+ 56.85 



2.50.3 

^ - 10.32 
3! 



-1.05^6 



a«(4.5) 



as (4.5) 



6! 



C - 24.16 



as (4.5) 

as(t) 



0.93 



as (4.5) 



0.58 ^2 
2[ 



-2.27 M 

^ + 3.918 
4! 



/ as(t) 
Us (4.5) 



-1.584^5 

5[ 



(5.19) 



where it is necessary to introduce a term cx a; in order to get a good description at high N. The 
full xP2L{x,t) is shown in fig. 16. a along with the two contributions above. The (5(1 — x) term is 
replaced in the power series by the x-dependence in the 0{as{Q^)) quark-gluon splitting function, 
i.e. x{x'^ + (l — x)^), normalized by 1.5 to give the correct AT — > limit. This corresponds to a slight 
modification of the usual physical anomalous dimension in terms of the 0{as{Q^)) longitudinal 
gluon coefficient function, but may be viewed as an analytic function with the correct N ^ limit 
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which aids presentation here.^° The d{l — x) terms in the power-suppressed contribution are very 
small, and are simply left out. In fig. 16.b we see xP2L{x,t) plotted as a function of x along with 
the naive LO BFKL calculation with coupling asiQ"^), and in order to illustrate the contribution 
of the higher-order terms, also the zeroth order contribution 1.5x(x^ + (1 — ^Y)- As with Pll{x, t) 
one can see that P2L{x,t) has a dip at small x before the eventual power-like growth sets in, again 
only for x < 10^^, and as with all calculated quantities the running coupling corrections severely 
diminish the strength of the small x growth. 

We can also try to investigate the effect of NLO corrections on physical quantities. In terms 
of partons the only known NLO correction is that to the gluon splitting function - there is simply 
no information on the NLO corrections to coefficient functions or the quark splitting functions. In 
terms of the physical anomalous dimensions, similarly there is no real information for r2L{N,t), 
but the situation is better for rLL{N,t). Let us look at the expression in terms of the partonic 
quantities (5.14), for the moment in the leading ln(l/a;) expansion without resummed running 
coupling corrections. At LO in 1/iV, r££(Ar, t) is equal to ^g^{N,t) since the differentiation of 
the log of the coefficient function with respect to t automatically introduces an extra factor of 
PoasiQ"^)- At NLO in 1/N r|ff'^(iV,i) picks up a contribution from ^^g^°{N,t) which is (largely) 
independent of the running coupling, and the contribution from the derivative of the LO coefficient 
function, which is entirely running coupling dependent. Hence, by knowing 7^^'^(7V, i) we know 
the whole of Tl^'^{N, t) before resuming running coupling corrections. Hence, we might hope that 
using an expression of the form (5.4), but corrected in the way described in the previous section 
for the NLO corrections to the kernel, we might calculate the full NLO, running coupling corrected 
BFKL expression for T ll{N ,t). Unfortunately, this is not quite the case. This can be appreciated 
by again using (5.14). When solving this NLO corrected expression for J'e,l{N ,t) one includes all 
the running coupling corrections to ^^J"'^ {N just by the manner of solving the equation. But 
without knowing the NLO correction to the coefficient function one misses a whole series of terms of 
the form as{Q'^){k)as{Q'^)T f{as{Q'^)/N) which would come form the '^'"^'^^jf^^'*^^ term.^^ Thus, 
we do not yet know the full running coupling corrections to the NLO contribution to Tll{N, t). 

I will proceed to calculate the "NLO" corrected TLL{N,t) on the assumption that since the 
resummation of the running coupling corrections stabilizes the perturbative expansion the missing 

This modification to the physical spUtting function will be discussed in a future paper. 
Some of these are automatically generated by using the NLO kernel in our solution, but the full set 
requires also the NLO correction to the hard scattering cross-section which will lead to NLO corrections to 

hL,g(,7)- 
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running coupling corrections will not lead to anything other than minor corrections. It is straight- 
forward to generalize the results of section 4 to the case of the physical quantity. Essentially we 
just replace (4.21) by 

(5.20) 

where we are currently missing a further term of the form 



Using (5.20) we can calculate both the power-series and power-suppressed NLO contributions to 
^LL{N,t) and hence PLL{x,t). The LO -|-"NLO" values of the intercept for the asymptotic power- 
like behaviour are shown in fig. 4. These lie very slightly below the LO-I-NLO intercepts for the 
gluon, and hint at perhaps a more rapid convergence for the physical than for the gluon. 
However, we would expect the missing contributions to lower the intercept a little more. The 
"NLO" corrected anomalous dimension r^1~^^^'^{N,t) is shown as a function of iV for f = 6 in 
fig. 14.b. It is very similar to that at LO until very low N where the difference in the leading 
intercept starts to become apparent. 

As for the NLO correction to xPgg{x,t) the power-series is not very convergent an to work all 
the way down to = IGeV^ and x = 0.00001 we again need the first 20 or so terms. Hence the 
power-series contribution is 

19 r rirna x / ^jji nn — ui — l \ 

xP™(«,(Q2),^) =«^(q2)^ E ^':iQ')[Knm^^, + K^sm^-x)], (5.22) 

n=l m=0 ^ ■ 

where the coefficients are listed in table 2. The power-suppressed contribution is 
(-0.183 exp(-0.51t)(5(l - x) + exp{-t) ( 31.90 ( i. I -80.22( "'^ ^ 



0:^(4.5)7 ■ V"s(4-5) 

+ 56.67( ^)°"|+S.OI7(^)"'|-25.925(^)"-"'| (5.23, 
\a,{4.5)J 2! V««(4-5)y 3! V««(4-5)y 4! 

+ 10.28I^V"t-1.298/'"-<"^"°«° 



as(4.5) J 5! ■ V««(4-5) / 6 

The NLO correction to the splitting function xP^l'^{x,t) is shown, minus the contributions oc 
(5(1 — x), in fig. 17.a. Clearly there is a very large cancellation between the power-series and power- 
suppressed contributions resulting in a relatively small total NLO correction. We can see that 
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unlike for the the gluon this NLO correction is actually positive in some regions of rather than 
everywhere negative. We also see from fig. 17.b that the NLO splitting function is quite similar to 
the LO splitting function over the whole x range. 

However, as with the gluon, the real test of perturbative stability is the evolution of the 
structure function itself. This is shown in fig. 18 where the evolution of a model for the structure 
function Fl{x, Q^), i.e. (1 — is shown both for the LO running coupling splitting function, 

and for the "NLO" corrected one (all 5{1 — x) contributions other that at first order in as{Q^) are 
included). Also shown is the evolution due just to the double-leading-log term P{x) = as{Q^)/x. 
Compared to the evolution of the gluon shown in the previous section we see that the additional 
running coupling contributions due to the t-derivative of the coefficient function have slowed the 
LO evolution below that of the double-leading-log result over the whole range of x (except very 
high x), and this will only cease to be true at very small x indeed, when the power- like growth 
of the physical splitting function finally sets in. In this case, however, the difference between 
LO and LO-|-"NLO" is much smaller than for the gluon, and the perturbative expansion seems 
very stable indeed. As with the NLO corrections to the intercepts this might be a sign that the 
expansion converges more quickly for the physical structure functions than for the unphysical gluon 
structure function. However, as a note of caution, the missing contributions at NLO are likely to be 
negative in general, and this difference between LO and NLO evolution will probably be increased 
a little. In fact it is desirable for these missing contributions to be non-negligible. Whilst if we 
decrease t to 4.5, i.e. IGcV^, at NLO everything remains relatively stable for the gluon, the 

physical splitting function starts to develop extreme behaviour at this low scale - the minimum at 
x ~ 0.01 becomes much lower and the peak at x ~ 0.0001 becomes very much higher. This trend 
is illustrated in fig. 11 b. which shows the rcnormalization scale dependence of P^^^''^^'^ {x,t) 
for t = 6. Clearly there is very good stability for an increase in scale, but it is not so good for a 
decrease in scale (though since the splitting function oscillates, the variation washes out to a large 
extent when evolution is calculated). There is very good stability in both directions if one examines 
the variation for a slightly higher t, say t = 8 {Q^ ~ 30GcV^). This instability in the physical 
splitting function results in instabilities in the evolution at t = 4.5, even though it appeared to be 
very stable at t = 6. Hopefully, the inclusion of the missing terms will help stabilize this evolution, 
though it may simply be a sign that at this low some non-perturbative contribution is becoming 
essential. 

6. Conclusions. 

In this paper I have shown that it is possible to obtain analytic solutions to the LO running 
coupling BFKL equation for the Q^-dependent parts of the gluon structure function and for the real 
physical structure functions F2{x, Q"^) and Fl{x,Q^). This results in a resummation of the leading 
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ln{l/x) terms at each power in as{Q^) and also of the leading powers in Pq at each power of as{Q^) 
and ln(l/x). However, the Qg-dependent gluon input is plagued by contamination from infrared 
nonperturbative physics, and has an inherent ambiguity of 0(A^/Qq). The analytic expressions 
may be expressed in form of a power-series in as{Q^)- In practice the main features of the solution 
are almost completely determined by only the first handful (~ 5) of terms in the expansion, in 
complete contrast with the case of fixed coupling, where an all orders summation is needed. In fact 
the perturbative series for the structure functions is not convergent, and the analytic expression 
is most accurately obtained by this truncation. The small remainder, which roughly speaking 
is suppressed by powers of (A^/Q^), may be calculated from the difference between a numerical 
solution with the analytic solution, and then modelled by an analytic expression of and N, 
which may easily be transformed to x-space. There are two points to note here. Firstly, this power- 
suppressed condition is both well-defined and is nothing to do with higher twist operators. Even 
though there are infrared (and ultraviolet) renormalons in the untruncated perturbative expansion, 
they only appear due to the impossibility of expressing the Q^-dependent part of the the structure 
functions as a power series in as{Q^), not because of some inherent ambiguity at leading twist 
as is often the case with renormalons. Hence, they are circumvented completely by this manner 
of calculation. Secondly, this procedure of an analytic calculation as a truncated power series 
plus a numerical calculation of the power-suppressed part, which is then modelled, seems to allow 
for the most accurate determination of x-space quantities. Transformation of numerical moment 
space expressions to x-space are subject to errors, and the magnification of the power-suppressed 
contributions in x-space, compared to moment space, seen in this paper highlights the potential 
effect of small errors in moment space when ultimately working in x-space. Hence, obtaining as 
accurate an analytic moment space expression as possible is vital in ultimately obtaining good 
accuracy for splitting functions and the evolution of structure functions. 

It is also demonstrated that there are well-defined, calculable higher-twist contributions due 
to the transverse degrees of freedom of the two-gluon operator. However, both the normalization 
and splitting functions of these genuinely higher twist operators decrease quickly as x ^ (roughly 
like x°'^ cos(0.51n(l/x))) when the small x resummation is performed. Unlike leading twist, this 
is largely insensitive to the running coupling corrections. This result is only apparent from resum- 
mation, and a fixed (small) order in as{Q^), particularly first order only, gives very misleading 
results. Hence, this one form of higher twist does not lead to any sizable correction at all at small 
X and Q^. It is possible that this unambiguous, small-x vanishing higher-twist contribution to the 
two-gluon operator is responsible for the absence of a genuine ambiguity in the leading twist anoma- 
lous dimensions. However, I note that this paper has nothing to say about the size of shadowing 
corrections coming from four gluon operators, except to point out that the double-leading-log type 
calculations often performed are likely to lead to huge overestimations. Neither does it consider the 
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power-suppressed corrections due to nonperturbative effects wliicfi mix with higfier twist, leading 
to mixing with leading twist, and may well be important at extremely small x [15] [29]. 

The calculated expressions for leading twist structure functions may be used to produce LO 
expressions for the splitting functions and coefficient functions for physical processes, and also the 
physical splitting functions which allow one to work directly in terms of physical quantities. My 
results prove that the effect of the running of the coupling is to weaken the asymptotic power- 
like growth of the splitting functions severely compared to the naive BFKL results, and even to 
lower the splitting function below the as{Q^)/x contribution for 0.001 ^ x ^ 0.2. It is also noted 
that the asymptotic behaviour of the form is often not approached even approximately until 
X « 0.00001, with the required x decreasing with increasing Q^, and is therefore by no means a 
good indicator of physics at present or future colliders. In fact it is very likely that unitarization 
will stop this true power-like behaviour ever being seen. Rather than the intercept, the detailed 
expressions for the splitting functions and coefficient functions are needed in order to really calculate 
the evolution at realistic values of x. 

The procedure can also be extended to NLO without any real modification, though there is 
some ambiguity in precisely what the best definition of NLO is.^^ The choice is made so that 
the expressions for the structure functions are genuinely only a single power of ag^Q"^) down on 
LO, up to (3oas{Q'^) corrections, but in j{N,t) = dlia{G{N,t))/dt the full NLO expression for 
{G{N, t))~^ is used, rather than truncating its expansion at NLO, and hence the full NLO correction 
to the intercept is obtained. This has little effect until extremely small x. Unlike leading ln(l/a;) 
calculations without resummation of running coupling effects the NLO correction to the gluon 
splitting function here is moderate, both for the value of the intercept and for the exact size of the 
splitting function and the evolution of gluon structure function for x > 10~^. Hence, this running 
coupling resummation does a great deal to stabilize the perturbative series. Unfortunately it is 
not yet possible to calculate the complete NLO correction to any real physical quantity, though 
one may come close for Pll{x, t), the splitting function governing the evolution of the longitudinal 
structure function in terms of itself, which is very similar to Pgg(x, t). In this case only a subset of 
the running coupling corrections to the NLO in ln(l/x) part are still unknown. For Fl the stability 
of the perturbative series looks even better than for the gluon as long as ^ 4GeV^, but begins 
to deteriorate below this, perhaps due to the missing corrections. 

Let me also comment briefly on other methods which attempt to incorporate the NLO correc- 
tions (and beyond) to the BFKL equation. Firstly I note that my previous conjecture that the effect 
of the running coupling in the BFKL equation could be accounted for using an x-dependent scale 
for the coupling [20], resulting in falling coupling for decreasing x, turns out to be essentially correct 

The power-series expressions also become very complicated at NLO. It will probably ultimately be 
more convenient to model them accurately with some simpler function of x and t similar to the manner in 
which the power-suppressed contributions are treated at present. 
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so long as the change in the scale of the coupling is moderate compared to the scale itself, though 
it fails if this condition is not satisfied. In practice this condition is identical to that specifying that 
diffusion in the fixed coupling BFKL equation is not too large, and therefore that the virtualities 
sampled in the running coupling equation are not too far away from . This results in the require- 
ment that ^ 201n(l/x) [35]. This is true for all but the lowest x and at HERA. I also note 
that my approach is completely consistent with that in [14] [15], with both being built upon the 
running coupling BFKL equation essentially introduced long ago [23] [24] [17] [18] and generalized 
beyond LO in [14]. The differences to this approach are that I ignore the coUinear resummation 
which is a central theme in this work, since as I stress it is an unnecessary complication in the 
calculation of splitting functions - the running coupling effects being the most important and stabi- 
lizing the calculation themselves; that I concentrate on solving very accurately and precisely for the 
Q^-dependent part of the gluon and structure functions, obtaining splitting functions over the range 
of X and relevant for a phenomenological treatment; and that I also ignore the complication of 
a real regularization of the coupling in the infrared region (this latter point is also considered in 
[36]). Hence, I obtain detailed accurate results for all splitting functions and coefficient functions 
in closed form, but ignore contributions considered in these papers which are necessary if investi- 
gating single-scale processes and/or potential nonperturbative effects (which may be important for 
splitting functions at low and very small x [29]). There is less similarity with other approaches. 
Even though that in [37] claims to in some sense be dealing with the scale appropriate for the 
coupling in this problem, it has no overlap with the approach in this paper, and comments on this 
approach can be found in [20] . Also there is no connection with the approach in [38] which adopts a 
phenomenological approach to resummation beyond fixed orders in ln(l/a:) in terms of the asymp- 
totic powcr-likc behaviour, which is a free parameter, and which consequently loses true predictive 
power for the evolution at small x. Finally, there also seems to be no overlap with the approach in 
the first part of [39] which incorporates subleading effects via a kinematic constraint while solving 
the BFKL equation - resulting in an anomalous dimension which includes a resummation of some 
subset of higher order contributions, none of which is concerned with the running of the coupling, 
but which stabilizes the calculation. (The latter part of [39] also includes a running coupling and 
infrared regularization, but concentrates on the normalization rather than the evolution.) In this 
sense it has some similarities to the resummation of coUinear logs in [13], which also stabilizes 
results even with fixed coupling (and which is essential in single scale processes). Hence, there 
appear to be a number of ways in which the apparent poor convergence of the perturbative series 
at small x can be improved. However, since one must ultimately deal with the contribution of the 
running coupling in all perturbative QCD calculations I prefer to concentrate on this feature and 
consider just the resulting resummation combined with the ln(l/a:) resummation, which results 
in explicit results in terms of an ordered power-series in the well-defined quantities ^^(Q^), ln(l/x) 
and /?o- This stabilizes the small x expansion without consideration of these other effects, indeed 
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it leads to the most divergent terms as a; — > [20] and alters the complete singularity structure, 
and moreover is easy to directly incorporate into the usual calculation of partons and structure 
functions in terms of the coefficient functions and splitting functions. 

It will, of course, be interesting to examine the effect of incorporating my resummed corrections 
to splitting function in a global fit to structure function and related data. Such an analysis will 
also need to include a precise explanation of how the small x relevant expansions derived in this 
paper must be combined with the normal order-by-order in as(Q^) expansion, and potentially large 
ln(l — x) expansions. Full details of such a fit, and the complete procedure used, will appear in a 
future paper which awaits the release of new data from a number of experimental collaborations. 
Prom the analysis of presently published data it is clear that the quality of such a fit is improved 
by inclusion of these small x resummed corrections^^, and that the predicted F/,(,t, Q^) is smaller 
than that form a NL0-in-Q!s((5^) fit, but much more regular in shape at low than that seen in 
[5].^^^ This can be seen as a solution to the lack of convergence of Fl{x,Q^) apparent as one goes 
from LO to NLO to NLO in the conventional expansion scheme which is seen in [5]. 

Hence, I conclude by claiming that this paper outlines a method for including the most complete 
resummation of splitting functions (and coefficient functions) which is needed at small x, and 
satisfies the theoretical requirements of stability of the perturbative expansion and the minimum 
of model dependence as well as the more practical considerations of being in a closed form which is 
easy to implement. It will prove useful in an analysis of structure function data, and in a prediction 
of related quantities relevant for the Tcvatron and the LHC. However, at present it only really exists 
at LO (and not even that for many quantities), and for full implementation the calculation of the 
NLO impact factors within the BFKL framework is urgently needed. Once this is done, a truly full 
NLO analysis of structure functions, which will be equally valid over the full perturbative range 
will be possible. 
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A prediction for dF2{x, Q^) / din for preliminary ZEUS data is shown in fig. 35 of [40] which shows 
the clear improvement compared to a convention NLO-in-as((5^) treatment. 
Figures showing this can be seen in [41]. 
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Table 1. 

The coefficients Knm in 

xPj^^O(e,a.(Q2)) = a«(g2) J2 ^s{Q')[Kn J 2\ +Kn5m^-x)]. (6.1) 

n=l m=l ^ ' 



n 


m = 5 


m = 4 


m = 3 


m = 2 


m = 1 


m = 




m = 11 


m = 10 


m = 9 


m = 8 


m = 7 


m = 6 




m = 17 


m = 16 


m = 15 


m = 14 


m = 13 


m = 12 




m = 23 


m = 22 


m = 21 


m = 20 


m = 19 


m = 18 


1 












-0.4236 


2 










-1.354 


1.611 


3 








-7.000 


30.22 


-34.63 


4 






-5.686 


46.92 


-103.2 


63.85 


5 




-16.14 


193.5 


-797.2 


1373 


-918.0 


6 


186.0 


-971.0 


2518 


-3323 


2045 


-458.9 














-14.35 


7 


-1386 


5051 


-9865 


10113 


-4281 


709.2 












-10.60 


192.0 


8 


21431 


-59800 


99225 


-95325 


49058 


-11483 










-24.48 


511.5 


-4497 


9 


70532 


-46099 


-25896 


59631 


-29684 


2798 






-17.21 


349.4 


-3100 


15284 


-44034 


10 


-126387 


-261087 


735693 


-761882 


373984 


-77690 




-12.01 


326.1 


-3758 


23801 


-88010 


179647 


11 


8688676 


-9665206 


6981022 


-3087487 


771318 


-102010 




1117.6 


-15044 


119789 


-620744 


2179220 


-5256680 














-37.57 


12 


-1.621x10^ 


1.864x10^ 


-1.288x10^ 


5044618 


-962638 


64963 




46536 


-211318 


563130 


-547416 


-1766356 


8225690 










-18.36 


506.4 


-6353 


13 


1.139x10^ 


- 1.236x10^ 


8.576x10^ 


-3.677x10^ 


9106782 


-1186015 




-400076 


1467162 


-2855626 


-645255 


2.136x10^ 


-6.623x10^ 








-12.58 


453.6 


-7298 


68149 


14 


4.425 xlO* 


-2.089x10^ 


6.225x10^ 


-1.129x10^ 


2061380 


-386008 




1.425x10^ 


-5.541x10^ 


+1.589x10* 


-3.345x10* 


5.135x10* 


-5.625x10* 






-46.85 


180.1 


-32842 


359398 


-2678626 


15 


-7.723x10^ 


1.340x10^ 


1.391x10^ 


-1.075x10^ 


3.170x10^ 


-4331143 




-1.146x10^ 


1.108x10^ 


-4.028x10^ 


9.069x10® 


-1.361x10'' 


1.339x10^ 




639.2 


-103283 


101509 


-640342 


2523533 


-4527424 














-18.33 
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n 


m = b 


m = 4 


m = 6 


m = 2 


777 = 1 


777 = 




171 = 11 


m = lU 


A 


m = 8 


777 = 7 


777 = 6 




Tfl = 17 


?Tl = 16 


m = 15 


771 = 14 


777 = 13 


™ 19 

777 — IZ 




m = 16 


m = iz 


m = zl 


777 = 20 


777 = 19 


777 = lo 


lb 


1.4U3X1U 


O A£?/^ V > 1 a9 

-8.066x10 


O 1 TO V , 1 a9 

3.172x10 


-8.188x10 


1.324x10 


-1.227x10 




O AO'? V , 1 a9 

z. 987x10 


1 /I A A V , 1 a9 

-1.490x10 


/I C AO V , 1 a9 

4.598x10 


A 0£? C V > 1 a9 

-9.865x10 


1 CO A V , 1 aIO 

1.529X 10 


1 TO /I V , 1 aIO 

-1.724x10 








1 n7f>Q/iQ 
-iu / Do4y 


O'iOi.O'iK} 


-i.4tol X lU 


-yui4tDy4i: 












-lz.3 / 


552. i 


17 


-1.108x10 


7.809x10 


-3.096x10 


6.807x10 


r-7 -| o > , 1 a7 

-7.184X 10 


1672556 




y.oObxlO 


1 \ / 1 nlO 

-1.676X 10 


O 1 'TO \ / 1 aIO 

z.l7zxlO 


-1.8z4x 10 


a A c \ / 1 a9 

6.057x10 


6.571 X 10 




ouoyuo 


-it OO40Z 


0.044X iU 


-z.yoz X iu 


i.zZ / X iU 


-o.ooy X iu 










-57.64 


2645 


-ooOoo 


18 


7.129x10 


-4.890x10 


2.107x10 


TO o > - 1 aQ 

-5.733x10" 


A r o 1 , - 1 aS 

9.531x10 


O O O yl , - 1 aT 

-8.884x10 




-3.690x10 


6.574X 10 


O AO /I ^ , 1 a1 1 

-8.024x10 


5.604X 10 


O 1 A CT ^ , 1 aQ 

3.195x10 


r TTo ^ , 1 a1 

-5.773x10 




1 uooz±y 


9 orifi V 1 


A 4fi7 V 1 Ci^ 


7 ^1^^ V 1 


zi 1 V 1 n9 

loo A lU 


1.40U A lU 






-17.59 


741.8 


-14788 


180754 


-1460081 


19 


4.467x10^1 


-1.573x1011 


3.986x101° 


-7.340x10^ 


9.737x10* 


-8.398x10^ 




4.962x10^1 


-9.477x1011 


1.402x1012 


-1.602x1012 


1.399x1012 


-9.198x1011 




-4.778x10'^ 


-6.248x10^ 


1.948x10^ 


-1.372x101° 


6.153x101° 


-2.005x1011 




-11.75 


620.5 


-15217 


225792 


-2200274 


1.398x10^ 



The scries for the part proportional to 6{1 — x) is more convergent in as{Q'^) and for all 
^ IGeV^ is given accurately by 



a,(Q')5(l - x) [%.{){Poas{Q')f + 139.5(;Soa,(Q'))' + 38.88(;aoa. + 964.2(^o««(Q'))^ 
+ 167.0(^oa«(Q^))^ + 5605(^0". ((3'))'°] ■ 

(6.2) 
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Table 2. 

The coefficients Knm in 

xPil^O(e,a.(Q2)) = a«(g2) J2 ^s{Q')[Kn J 2\ +Kn5m^-x)]. (6.3) 

n=l m=l ^ ' 



n 
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m = 1 
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m = 11 
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m = 13 
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1 
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-1.354 
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3 
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-29.49 


4 






-5.672 
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-222.22 
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5 




-15.84 
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6 


243.4 


-1444 
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5206 
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-17.45 


7 


-5265 
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8 
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-84630 
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-187932 
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-6.545 


448.6 


-5158 


9 
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-122925 
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-20.94 
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-5027 


31574 
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10 
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1019004 
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-22.36 


814.7 


-12094 


91607 
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11 


3.965x10^ 


-5.343x10'^ 


4.627x10^ 


-2.425x10^ 


7013331 
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1291 


-23513 
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-1492972 
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-1.956x10^ 














-20.00 


12 


-2.424x10^ 


2.970x10^ 


-1.702x10*^ 


1.214x10* 


-3.553x10^ 
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60402 


-153923 


-952951 


1.057x10^ 


-4.805x10^ 


1.328x10* 










-22.33 
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-8519 


13 


-1.514x10^ 


- 7.7974x10^^ 


1.993x10^ 


-1.426 xlO^ 


4.814x10^ 
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1.199x10^ 


-4.586x10^ 


1.252x10^ 


-2.356x10® 


2.796x10* 








-23.25 


1171 
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14 


4.934 xlO^ 


-2.902x10^ 
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-1.227x10* 


-1.883x10^ 
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5.595x10^ 


-2.675x10^ 
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-2.317x10^ 
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-33.40 


2462 


-61331 


876031 


-8361477 


15 


-2.576x10^° 


1.702x10^° 


-7.397x10^ 


1.968x10^ 


-2.961x10® 


2.201x10^ 




-3.263x10^ 


1.473 xlO*^ 


-4.775 xlO^ 


1.147x10^° 


-2.052x10^° 


2.708x10^° 




781.7 


-12818 


114646 


-395405 


-2857628 


4.706x10^ 














-22.34 
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The scries for the part proportional to 6{1 — x) is more convergent in as{Q'^) and for all 
(52 ^ lGeV2 is given accurately by 



a,(Q')(5(l - x) [-0.3094(^oa«(Q')) - 3.856(^oa«(Q'))' + 6.376(^oa,(Q'))^ - 50.36(;9o««(Q'))^ 
+ 340.0(;9oa.(Q'))^ + 55.51(^oa.(Q'))' - 1600(;9oa.(<5'))^ + 2838(^oa«(Q'))^ - 8457(;9oa«(Q'))^ 
+ 24526(;9oa,(Q2))i° + 57602(^oa.(Q'))" - 325984(^oa.(Q'))'' + 477536(^oas(Q'))''] • 

(6.4) 
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Figure Captions 



Fig. 1. The branch points and cuts associated with exp(— Xo(7)/(/3o^)) and the saddle-point 
contour, the Gamma-function contour and the numerical integration contour. 

Fig. 2. The expression A4(— l/(^o-^)) as a function of N compared to l/(^o-^)^- 

Fig. 3. The Q^-dependcnt part of the gluon structure function, QE{N,t), and of dQE{N,t)/dt 
as a function of for t = 6 {Q^ '-^ 6GeV^). The Q^-independent factor of 
— sm{TT / {j3oN))T{—l/ {l3oN)) exp{—^E / {l3oN)) is included in both in order to produce 
a smoother A^-dcpcndent normalization of the functions. 

Fig. 4. The positions of the leading poles in the anomalous dimensions for the gluon structure 
function at LO and NLO, and for at LO and NLO. 

Fig. 5. a. The anomalous dimension for the gluon structure function at LO plotted as a function 
of for t = 6 ~ 6GcV^). Also shown is the 0(as(Q^)) contribution as{Q'^)/N, and 
the full naive LO BFKL anomalous dimension, b. The anomalous dimensions for the 
gluon at LO and at NLO plotted as functions of N for t = 6. 

Fig. 6. a. The splitting function xPg^(x) and its power series and power-suppressed contributions 
plotted as functions of a; for t = 6. b. The splitting function xPg^{x) plotted as a function 
of X for t = 6 (Q^ ~ 6GeV^). Also shown is the 0{as{Q'^)) contribution ^^(Q^), and the 
naive LO BFKL splitting function with coupling as(Q^). 

Fig. 7. The next-to-leading twist splitting function for = 0.2. 

Fig. 8. The value of exp(— + Xq{'^) / {J3qN)) , along the real axis for = 0.4 and Iq = 2, along 
with 7 = 7o^^(io), 7o^^(i) for t > to. 

Fig. 9. a. The splitting function xP^g^^{x) and its power-series and power-suppressed contribu- 
tions plotted as functions of x for t = 6. b. The splitting functions xPg^^^^'^{x) plotted 
as a function of x for t = 6 (Q^ ~ 6GeV^). Also shown is the 0{as{Q'^)) contribution 
«s(Q^)) and the LO contribution xPg^{x). 

Fig. 10. The values of dG(x, Q2)/dlnQ^ for G{x,Q'^) = x-°-^{l - xf, due to the LO splitting 
function Pg^{x) and the LO-I-NLO splitting function PgJ-'~^^^'~'{x), plotted as functions of 
X for t = 6 (Q^ ~ 6GeV^). Also shown is the evolution due to the 0{as{Q'^)) contribution 
P(x) = a,{Q^)/x. 

Fig. 11. a. The renormalization scale variation of the LO-I-NLO splitting function Pg^^^^^ . 

Shown are the three choices of scale Q^, 0.5(5^ and 2Q^ for t = 6, i.e. ~ 6GeV^. b. 
The same for the LO-hNLO physical splitting function p^o+nlo 

Fig. 12. a. The full leading ln(l/.T) plus running coupling corrections coefficient function 
xCL,g{x,t) plotted as a function of x for t = 6 and Nf = 4. Also shown are the con- 
tributions from the power-series and the power-suppressed part. Note that the term 
oc 5{1 — x) in the power-series is replaced by the full 0{as{Q^)) contribution 6x^(1 — x), 
and the terms oc (5(1 — x) in the power-suppressed part are absent, b. xCf^(x,t) plotted 
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as a function of x for t = 6 and Nf = 4. Also shown is the coefficient function obtained 
from the naive LO BFKL calculation, and the contribution at 0{as{Q^)) alone. 

Fig. 13. a. The full leading ln(l/a;) plus running coupling corrections coefficient function xPqg{x, t) 
plotted as a function of x for t = 6 and Nf = 4. Also shown are the contributions from 
the power-series and the power-suppressed part. Note that the term oc (5(1 — x) in the 
power-series is replaced by the full 0{as{Q'^)) contribution 1.5x{x^ + (1 — x)'^}, and the 
terms oc 5{1 — x) in the power-suppressed part are absent, b. xP^^ {x,t) plotted as a 
function of x for t = Q and Nf = 4. Also shown is the coefficient function obtained from 
the naive LO BFKL calculation, and the contribution at 0(a<j((5^)) alone. 

Fig. 14. a. The anomalous dimensions for the gluon structure function at LO and for Tl {N, t) 
at LO plotted as functions of for t = 6. Also shown is the 0{asiQ^)) contribution 
common to each. b. The anomalous dimensions for J^L{N,t) at LO and "NLO" plotted 
as functions of A/" for t = 6. 

Fig. 15. a. The splitting functions xP^j^{x) and its power series and power-suppressed contribu- 
tions plotted as a function of a; for t = 6. b. The splitting function xP^2{x) plotted 
as a function of x for t = 6 (Q^ ~ 6GeV^). Also shown is the 0{as{Q'^)) contribution 
Q!s(Q^)) the gluon splitting function Pg^{x) and the naive LO BFKL splitting function 
with coupling agiCp)- 

Fig. 16. a. The full leading ln(l/.T) plus running coupling corrections physical splitting function 
xP2Lix,t) plotted as a function of x for t = 6. Also shown are the contributions from 
the power-series and the power-suppressed part. In the power-series the part oc 6(1 — x) 
is replaced by 1.5x{x^ -|- (1 — x)^) while in the power-suppressed part this contribution 
is simply absent, b. The physical splitting function xP2L{x,t) plotted as a function of 
X for t = 6 along with the physical splitting function obtained from the naive LO BFKL 
calculation with coupling as{Q^) and the zeroth order contribution. 

Fig. 17. a. The splitting functions xP^l'^{x) and its power series and power-suppressed contribu- 
tions plotted as functions of x for t = 6. b. The splitting function xP^^^^^'~' {x) plotted 
as a function of x for t = 6 (Q^ ~ 6GeV^). Also shown is the 0{q.s{Q'^)) contribution 
Q.s{Q^)-, and the LO contribution xP^^{x). 

Fig. 18. The values of dFL{x, Q'^)/d\n Q^, for Fl{x, Q^) = a;-°-2(l - x)^ due to the LO splitting 
functions Pff (x) and the LO-^NLO splitting function Pff plotted as a func- 
tion of X for t = 6 (Q^ ~ 6GeV^). Also shown is the evolution due to the 0{as{Q^)) 
contribution P{x) = as{Q'^)/x. 
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